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A unifying approacli to competing quantum orders in generalized two-leg spin ladders is presented. Hidden 
relationship and quantum phase transitions among the competing orders are thoroughly discussed by means 
of a low-energy field theory starting from an SU(4) quantum multicritical point. Our approach reveals that 
the system has a relatively simple phase structure in spite of its complicated interactions. On top of the U(l)- 
symmetry which is known from previous studies to mixes up antiferromagnetic order parameter with that of the 
p-type nematic, we find an emergent U(l)-symmetry which mixes order parameters dual to the above. On the 
basis of the field-theoretical- and variational analysis, we give a qualitative picture for the global structure of the 
phase diagram. Interesting connection to other models (e.g. bosonic t-J model) is also discussed. 
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I. INTRODUCTION 



In the past two decades, quantum magnetism has been serv- 
ing not only as effective theories describing insulating phases 
of strongly-correlated electron systems but also as theoretical 
laboratories to look for and test new concepts. The discovery 
of high-temperature superconductors sparked the search for 
unconventional or exotic phases which are quite different from 
the ordinary ferromagnetic/antiferromagnetic phases. Despite 
the effort of many researchers in searching for novelty, it is 
by now well-known that, in two- or higher dimensions, an- 
tiferromagnetic phases are found provided that spin frustra- 
tion is not very strong'. In order to suppress antiferromag- 
netism and stabilize exotic phases^, various mechanisms have 
been proposed. One realistic example of such mechanisms 
would be multispin-exchange interactions. Such interactions 
are expected to be crucial for explaining unusual magnetic be- 
havior in ^He absorbed on graphite^. Moreover, it was re- 
ported that a certain amount of four-spin ring exchange would 
be necessary to account for neutron-scattering experiments 
for the parent compound of high-temperature superconductor** 
La2Cu04 and for a spin-ladder compound^ '' La6Ca8Cu2404i. 
Extensive numerical simulations caiTied out^ ** for the two- 
dimensional Heisenberg antiferromagnets with four-spin ring 
exchange found various phases with unconventional orders: a 
spin-nematic phase** and a spin-liquid phase with topological 
ordering'. 

On the other hand, various approaches have been pro- 
posed in electron systems to unify several (and sometimes 
quite different) competing orders^"'^ and succeeded in clar- 
ifying the nature of the quantum phase transitions among 
them'^ '"*. Usually, in those approaches, extended symmetries 
are adopted so that mutually competing order parameters may 
be transformed to each other A typical example would be the 
SO(5) theory^'"' for the competition between d-wave super- 
conductivity and antiferromagnetism, where the order param- 
eters of (i-wave superconductivity and those of antifeiTomag- 
netism are combined to form a unified order parameter quin- 
tet. For a one-dimensional geometry (two-leg ladders), it is 



known that even larger symmetries SO(8)'^ '^ and SO(6)''"'^ 
can emerge at low energies and be useful for the description 
of the electronic phases at half-filling and away from half- 
filling, respectively. In particular, the existence of an extended 
symmetry might provide a route to classify one-dimensional 
gapped phases^'^. 

Unfortunately, no systematic approach based on extended 
symmetries is known for unconventional phases found in spin 
ladders and it is desirable to construct such theories. As the 
first step along this line, we investigate here two-leg spin 
ladders with four-spin interactions, since they possess high 
enough symmetry to unify various competing orders. Without 
the four-spin interactions, the two-leg spin ladder has a finite 
spin gap in magnetic excitations and short-ranged magnetic 
correlations^'. Basic physics of the two-leg ladder can be un- 
derstood by considering a ground state consisting of almost 
localized dimer singlets on antiferromagnetic rung bonds (see 
Fig. 1) and low-energy excitations over it. For this reason, 
the spin-gap phase in the usual two-leg spin ladder sometimes 
is dubbed as rung-singlet- or rung-dimer phase. If we change 
the rung coupling to ferromagnetic, the rung bonds will be 
dominantly occupied not by singlets but by triplets. Then, 
the system is effectively equivalent to the spin-1 systems and 
the knowledges in those systems may apply^^"^''. When the 
four-spin interactions (say, ring-exchange) are switched on, 
the size of the spin gap will be reduced and finally at a cer- 
tain critical strength it even vanishes^^"^**. Large-scale numer- 
ical simulations^^'-'^ suggested that the model has a rather rich 
phase diagram. In particular, a spin liquid phase with scalar 
chirality ordering was found^^. Such a phase breaks both 
time-reversal and parity symmetries and has been discussed 
previously in the context of any on superconductivity^'"^^. A 
hallmark of the unconventional phases in the phase diagram of 
the two-leg spin ladder with ring-exchange interaction is that 
neither singlets nor triplets dominate over the others. Hence, 
the conventional approaches starting from the limit of strong 
rung couplings (whether ferromagnetic or antiferromagnetic) 
is not very convenient to explore the nature of novel phases 
stabilized by four-spin interactions. The main goal of this pa- 
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per is to fully describe the nature of the unconventional phases 
and quantum phase transitions among them by an approach 
based on an extended symmetry. 
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FIG. 1: (a): Two-leg spin ladder. The unit of constructing the model 
Hamiltonian (see section II) is a pair of spin-I/2s enclosed in a bro- 
ken line (called rung in the text). Any site can be labeled by the chain 
index (1,2) and the rung index r. (b) and (c): Two typical states ap- 
pearing in two-leg ladders. The usual rung-singlet phase (b) and the 
staggered-dimer phase (c) to be discussed in the text. Arrows denote 
spin-singlet pairs. 



The organization of the present paper reads as follows. In 
section II, we construct the ladder Hamiltonian by requir- 
ing rotational invariance and equivalence of two constituent 
chains. Our goal is to unify several unconventional phases sta- 
bilized by four-spin interactions. To this end, we shall use an 
enlarged symmetry SU(4) which contains the ordinary (spin) 
SU(2). We shall also pay particular attention to an interesting 
symmetry {spin-chimlity transformation) which is a special 
case of the above SU(4) and commutes with the spin rotation. 
This symmetry, which has been introduced^"'^-' in the context 
of the ring-exchange two-leg spin ladder, exchanges the Neel 
order parameter with the vector chirality (the order parame- 
ter of a p-type nematic^^"*). It will play a crucial role in our 
analysis and give an important clue to understand the global 
structure of the phase diagram. 

A low-energy field-theory analysis will be developed in sec- 
tion III. Although various field-theory approaches are known 
for two-leg ladder systems- most of them start from 
the limit of two decoupled chains and will not be suited to 
investigating the phase structure when four-spin interactions 
are by no means small. Instead of starting from two weakly- 
coupled = 1/2 chains, we shall take the SU(4)-invariant 
point, which is a special case of the lattice Hamiltonian, as the 
starting point. The rotational- and the spin-chirality symme- 
try are beautifully incorporated into our low-energy effective 
action. 

In section IV, the phase structure and unexpected high sym- 
metry among unconventional phases will be discussed with 
the help of one-loop renormaUzation group (RG) calculation. 
In particular, we shall find four dominant phases where SU(4) 
(SO(6), more precisely) symmetry is approximately restored 
in the low-energy limit. The quantum phase transitions among 
these dominant phases will be investigated in sections V and 
VI. The crucial role played by the spin-chirality symmetry in 
transitions among spin-singlet phases will be revealed. In this 
respect, a bosonization scheme based on U(l)xSU(3) sym- 



metry will be introduced to extract relevant low-energy de- 
grees of freedom which govern the phase transition. By using 
an effective theory for the low-energy fluctuations of a vector 
doublet composed of competing order parameters, we shall 
clarify how time-reversal symmetry is broken in one of these 
unconventional phases. 

Readers who are not interested in the detail of the field- 
theoretical analysis may skip sections III, IV, and VI. Our 
main results have been already published in part as a short 
communication . 

In order to supplement the field-theoretical analysis, we 
shall carry out a variational- and a strong-coupling analysis 
in section VII. As has been described above, the usual strong- 
coupling expansion starting from the limit of isolated rungs 
is not very convenient. Instead, we shall start from the limit 
of isolated plaquettes to successfully describe the competition 
among several quantum phases. On the basis of results ob- 
tained from these analyses, we map out the global phase dia- 
gram and discuss a connection to a spin-1 bosonic t-J model. 
This relationship might be useful for the possible experimen- 
tal realization of the exotic gapped phases stabilized by four- 
spin exchange interactions. Indeed, in standard two-leg lad- 
der compounds discovered so far, four-spin exchange interac- 
tions (10-20% of leg/rung interactions) are not strong enough 
to stabilize the spin liquid phase with scalar chirality order- 
ing in these systems. However, with the connection to spin- 
1 bosonic t-J model presented in section VII, one can ex- 
pect that this phase (together with other unconventional ones) 
might be feasible in ultra-cold bosonic atomic gases in opti- 
cal lattices^^. Finally, section VIII presents our concluding 
remarks and technical details will be presented in the Ap- 
pendixes. 



II. THE MODEL AND ITS SYMMETRIES 

A. Building Blocks of a Hamiltonian 

As it has been discussed in the introduction, enlarged sym- 
metries are very powerful in unifying various competing or- 
ders. Therefore, it is necessary first to identify the enlarged 
symmetry in our problem. Since we are considering spin- 
1/2 two-leg ladders, the maximal symmetry would be SU(4). 
The appearance of SU(4) in our problem is easily understood 
by noting that four states (| tt), I Ti>, I it), I ii» on 
a single rung (say, r-th rung in Fig. 1(a)) span the four- 
dimensional defining representation (4) of 811(4)'*''''*'. Let us 
denote S ~ 1/2 spins on the first- (upper) and the second 
(lower) chains by Si ^ and S2,r, respectively. Then, one of the 
standard choices of the 15 generators X; on the r-th rung is 

X'{r) = Sl^ , X'+\t) = SI, (z = 1, 2, 3) , 

X'{r)^2S'l,Sl, = Gab{r) (z = 7, . . . , 15; a, 6, = 1,2,3). 

(1) 

The explicit matrix expressions of these generators are given 
in Appendix A. 
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Si 




...^2 


-even 


Si 


r^S2,r 


...Z2 


-odd 


Si 


r X S2,r 


...Z2 


-odd 



Now we are at the point of constructing the Hamiltonian 
by requiring the invariance under (i) SU(2), (ii) Z2 (1 ^ 2: 
interchange of the two chains. See Fig. 1) and (iii) reflec- 
tion with respect to horizontal Unks {link parity). Since we 
are considering S = 1/2, we can safely restrict ourselves to 
second-order polynomials in Sf and 82- We can divide them 
into (a) a scalar, (b) (SO(3)-)vectors and (c) rank-2 symmet- 
ric tensors. Of course, the scalar is given by Si_r-S2,r and the 
followings are all vectors: 



(2) 



Although the last one seems an antisymmetric tensor, it be- 
haves like a (pseudo)vector in the spin space and serves as the 
order parameter of the p-type spin nematic^''*. On top of them, 
we have a symmetric tensor whose components are given es- 
sentially by Gab'^~' 

Qab{f) = '5"j.'52,r + ^l,r^2,r ~ 2 (^'''>('') ~^ G't,a(r)) , 

(3) 

which reduces to a set of order parameters of the 71-type spin 
nematic^^ when the dimer bonds are occupied by triplets. The 
point here is that all these operators are essentially the SU(4) 
generators (see Appendix A. 2 for more details). 

Now we proceed to constructing a Hamiltonian. If we con- 
sider only interactions involving two neighboring rungs (i.e. 
four spins), the SU(2)- and Z2-invariance strongly restrict the 
possible interactions and we are left with the following seven 
ones; 

Ti-l = ^(Sl,r + S2,r)-(Sl,r + l + S2,r+l) (4a) 

r 

'H2 ~ ^(Si^r — S2,r)-(Sl,r+l — S2,r+l) (4b) 



^3 = 2^ I] Qab{r)Qab{r + l) 

r a,b=l 

= 4 [(Si^r'Sl,r+l)(S2,r'S2,r+l) 



(4c) 



+ (Sl,r'S2,r+l)(S2,r'Sl,r+l)] 

Hi ~ ^4(Si,r X S2,r)-(Sl,r+l X S2,r+l) 

r 

= 4 [(Si.T--Si^.r+l)(S2,r'S2,r+l) 

T 

— (Sl,r'S2,r+l)(S2,r'Sl,r+l)] 
7^5 = 2 [(Si.r — S2,r)-(Sl.r+l X S2,r + l) 



(4d) 



(4e) 



+ (Sl,r X S2,r)-(Sl,r+l — S2,r+l)] 
'^6 = ^^(Sl,r-S2,r + Si_r+1 •S2,r+1 ) (4f) 



'H'J — y~^(Si^r'S2,r)(Sl,r+l'S2,r+l) 



(4g) 



where the summation ^ ^ is taken over all rungs of the ladder. 
Aside from the four-spin terms 7^3, 7^4, and 7^7, we have a 
three-spin term H.^ which explicitly breaks time-reversal sym- 
metry. If S" = 1/2 comes from the electron spin (more gener- 
ally, magnetic moment of charged particles), three-spin term 
7^5 may result from the electron hopping on each plaquette. 

It would be useful to rewrite Hi and H3 in terms of spin-1 
operator defined on the r-th rungs; 



Hi — Tr'Tr+l , 

r 

^3 = E 



2 3 
Tr-Tr+1 + 2 (Tr-Tr+l) — iHe - - 



(5) 



In the above equations, the projection operators onto the 
triplet subspace Ptripiet(7') = Si^r-S2,r + 3/4 are impUed, i.e., 

Tr = Ptriplet(f)(Sl,r + S2,r )^'triplet (t") ■ 

These two blocks describe the interaction between effective 
spin-1 objects, that is, they dictate the magnetic part of the 
Hamiltonian. All the above seven interactions are used to con- 
struct the following general Hamiltonian; 

H = AHi+BH2+CH3+DHi+EH^+FH&+GH7 . (6) 

Of course, another set of interactions could have been used. 
For example, basis Hi ± 7^2 and H^ ± Hi could have been 
chosen instead of Hi, H2, H3, and Hi- In fact, the lat- 
ter choice is convenient when discussing the systems with 
SU(2)xSU(2)-symmetry (e.g. the spin-orbital model'*'''*^^-'*'*), 
while our choice here is useful when dealing with the spin- 
chirality transformation which will be introduced in the next 
subsection. 



B. Spin-chirality Transformation 

1. Construction 

As mentioned in the previous sections, the largest symme- 
try of the problem is SU(4), whose defining representation 
is spanned by the four states (singlet and triplet) on a single 
rung. As a subgroup of SU(4), there is an interesting U(l) 
symmetry^^ called the spin-chirality transformation. 

Let us consider two spins (Si and S2,r) on the r-th rung 
and look for a local unitary transformation Ur (0) which com- 
mutes with the SU(2)-rotation generated by Si^^ + S2,r- The 
commutation relations 



Si, 



»2,r 







(7) 



satisfied by the generator of Ur{0) uniquely (up to a constant 
phase) determine the following form; 



Ur{0) ^ CXp [i9 Ptriplet(?')] 



exp 
1 



— + Si_j.'S2,r 



(8) 



(l + 3e*'') + (e*''-l)Si,,-S2,r 
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In the above expression, we have adopted a sUghtly different 
definition from the original one in Ref. 33 for reasons which 
will become clear later By construction, it is obvious that 
Ur{0) is the only U(l) transformation which commutes with 
the spin-rotation symmetry. The U(l) transformation 
has remarkable properties. By fully utilizing the properties of 
S = 1/2, we can show that the following equations hold: 



algebra, we obtain the following rules: 



(total-spin conservation) 

Ur{e){Si^r-S2,r)Ul{e) 

= (Si.r — S2,r) C0b9 — 2(Si,r xS2,r) sillfl 

~ -(Sl,r — S2,r) Sn\9 + (Si_r X ^2,r) C0s9 



(9a) 



(9b) 



The first line (9a) is a trivial consequence of the requirement 
(7). The local U(l)-transformationZ//r(^^) can be readily gen- 
eralized to the whole lattice: 



(10) 



and all the above properties are preserved for L{{9) as well. 
In what follows, we shall call U{9) spin-chirality transfor- 
mation, since, as can be seen in the above equations (9b), it 
mixes up the antiferromagnetic order parameters Si — S2 and 
the vector chirality (or the order parameter of the p-type spin 
nematic-''*) Si x S2- 

For our purpose, it is helpful to view U {9) as an SU(4) 
transformation rather than as a non-linear transformation for 
two spin operators Si and S2. The latter two lines (9b) sug- 
gest that two antisymmetric (in 1 2) quantities Si — S2 and 
2(Si X S2) behave as an 0(2)-doublet under the spin-chirality 
transformation Z// (6'): 



Si — S2 
2Si X S2 



/ cos6' — sin6' 
I sin6' cos9 



Si — S2 
2Si X S2 



(11) 



The remaining 9 generators {Si 
under (6*). 



S2 and Gab) are invariant 



2. Duality 



Now we are at the position to discuss on the effect of spin- 
chirality transformation on our building blocks. After some 



n2^\{n2 + 'Hi) 



^(^2 + 7^4) 



(12a) 



1 



■ COS 29 {J-L2 -n^)-- sin26' Ti^ 



(12b) 
(12c) 



1 



1 



cos 261 (7^2 - + 7: sin26l T^g 



^ sin26l(7i:2 - W4) + cofi29n^ 



(12d) 
(12e) 
(12f) 
(12g) 



Now the reason why we have decomposed the Hamiltonian 
into Til, ■ • ■ can be easily understood from the above 
equations. The spin-chirality transformation for 6* = ^ is par- 
ticularly simple: 

7^2 ^ 'Hi , 7^5 I— > —7^5 (all the others are invariant) . 

(13) 

Hereafter, this special case 



V=U{9^ 7r/2) 



(14) 



will be called duality transformation^^'^^ , although it is rather 
different from the standard 'duality' which maps local objects 
onto non-local ones and vice versa. It readily follows from 
Eqs. (12a-12g) and (13) that a model with B = D and S = 
is invariant (or self -dual) not only for 6* = 7r/2 but also for 
arbitrary values of 9^^. In what follows, we will see that this 
enhanced U(l)-symmetry aX B = D will play a crucial role. 



C. Special Cases 

Before discussing the duality property of them, we mention 
interesting special cases of Hamiltonian (6). 

Ordinary two-leg ladder: Much is known^''-'^'^'^ about the 
ordinary two-leg ladder defined by the following Hamiltonian: 

H = J (Si^r'Sl^r + l + S2,r-S2^r+l) + J± Si^r'S2,r 

r r 

(15) 



The basic picture of the ground state is provided by putting 
dimer singlets on rung (i.e. J±) bonds and low-lying excita- 
tions may be understood as propagating dimer triplets"*^. 

Spin-orbital model: The Hamiltonian of the spin-orbital 
model is defined by"*' '*^ "*'* 

"Hso ~ J (Si^r'Sl,r+l + S2^r'S2,r+l) 



K (Si,r'Sl,r + l) (S2,r'S2,r+l) 



(16) 
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and is obtained by choosing 

2 8 

It was shown by weak-coupling analysis^^ and later by ex- 
plicitly constructing the exact ground state^^ that the model 
displays a staggered dimer (SD) ordering in a certain region 
{0 < K < 4 J) of the parameter space. 
If we further impose the restriction 



A = B ^ C = D { = 



the model reduces to the so-called SU(4)-spin-orbital 
model*^'*^''^^ . In this particular case, the Hamiltonian can be 
written as 



rGmng 



r^rung ^ 
15 



(17) 



where P(S,T) denotes a permutation operator for two S = 
1/2 modules (corresponding to S and T). This is an SU(4) 
generalization of the S = 1/2 Heisenberg model since 
P(Si,r, Si,r+i)-P(S2,r, S2,r+i) is nothing but the SU(4)- 
permutation operator if we regard the four states on a rung 
as 4. The model 7isu(4) is integrable"^^ and will be used as a 
starting point of the following analysis. 

As shown in Ref. 48, a term GHq can be added to 7isu(4) 
without spoiling integrability. According to the Bethe ansatz 
results'*^, we have two critical values of G; when G < Gc.i = 
7r/(2V3) + In 3/2, all rungs are occupied by triplets and the 
system is described by two gapless bosons while for G > 
Gc,2 = 4 the system is in the so-called rung-singlet phase 
and all excitations are gapped. 

Self-dual models: An important class of models is defined 
by the following choice of parameters: 

B = D, E = 0, A,G,F,G ^ arbitrary . 

This defines a family of models which are invariant under the 
full spin-chirality rotation U (6). Hereafter, we call this family 
of models self-dual models and the manifold characterized by 
the above set of parameters self-dual manifold. Obviously, the 
self-dual models have SU(2)spin x U( 1 )spii,.(.j,ii-ai symmetry. 
SU(3)xU(l )-models: On a special sub-manifold 



C 



1 



-Ji,B = D = -J2,E = 0,E,G = arbitrary 



of the self-dual manifold (obtained by setting yl = C in the 
self-dual models), the spin-SU(2) gets enlarged to SU(3) (see 



Appendix A 2). The conditions A = G and B = D are crucial 
for the SU(3)-invariance. In fact, these SU(3) and U(l) are 
broken simultaneously (SU(3) SU(2), U(l) t-^ Z2) when 
we move away from the self-dual manifold (B = D). 

composite-spin model: The so-called composite-spin 
model is defined by 



49,50 



"^composite " [(^l.r + ^2,r) ' (Sl,r+1 + S2,r+l) 

r 

-/9[(Si,, + S2,.)-(Si^,+i + S2,,-+i)]'] 



(18) 



which can be rewritten as 

'^composite = ^1 + -/3^ Hi — ^PH^ ~ S/JHg • 

This preserves the total spin Si.,. + 82,7- on each rung and 
the problem reduces essentially to that of a collection of finite 
chain segments. 

Note that this is a special case {B ~ D ^ Q) of the self- 
dual models. If we choose (3 = —1, we get A ^ G = 1/2 
and as a consequence we have SU(3) symmetry. This is in 
agreement with the well-known fact that the /3 = — 1 bilinear- 
biquadratic chain is SU(3)-invariant (the SU(3) Uimin-Lai- 
Sutherland model"*^'^'). 

spin-ladder with four-body cyclic exchange: The four- 
body cyclic exchange on elementary plaquettes made up of 
two rungs r and ?- + 1 can be recasted as 



Wcyc - (^'i^^' '•+!)+ Pi^i^^ ^ + 1)) 



Hi + + 2^6 + 471:7 + const. 



(19) 



where Pi{r,r + 1) and P^'^{r^r +1) respectively make a 
cyclic permutation and its inverse on a plaquette formed by 
rungs r and r + 1. Using this, the Hamiltonian for a two-leg 
ladder with four-spin exchange is given as^-^^: 

Hladder+4-spin 

r+l + S2,r-S2,r+l) + K^Hcyc 



(20) 



= ( ^ J + /V4 ) Hi + i J7^2 + KiHi 



+ ( Jr + 2/Q)H6 + 4X4^^7 ■ 

Note that the model with = J/2 is self-dual (B = D) and 
we can find not only the exact rung-singlet ground state but 
also the exact one-magnon state^'' (in notations used in Ref. 
26, Jiing = 2K4) for certain choices of parameters. The phase 
diagram for Jr = 1 has been mapped out in Refs. 26,27,29, 
30. 

We summarize the relation between parameters and the 
symmetries of Hamiltonian (6) in Fig. 2. First of all, in most 
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cases perturbation around SU(4)-point explicitly breaks the 
SU(4) symmetry. The fate of the system after the breaking of 
SU(4) is different according to whether or not the system is 
invariant under the spin-chirality rotation U(9); if the system 
is invariant (i.e. B = D), then we can have a high symme- 
try like SU(3)xU(l) (when A = C) or SU(2)xU(l) (when 
A 7^ C). These cases will be treated in section V. Otherwise, 
the system generically assumes the lowest possible symmetry 
SU(2) X Z2. 

i-orbital mode! A=B, C=D 



SU(2)XSU(2)XZ2 




SU(2)XZ2 











^elf-dJ^l 


C;tB-D 
► 


SU(3)XU(1) 


SU(2)XU(1)| 


F^O 





r ^A-C;tB = D| 1 A:+r' I 1 

I SU(4) I ► SU(3)XU(1) 1-k SU(2)XU(1)I 

^^___loi' 1"^^^^' ' ' ^ ' 

\ |B<D/^fZ):B^D 

.xsu(2)xz2 4 SU(2)XZ2 



SU(2)XSU(2)XZ2 


4 


SU(2)XZ2 


^=D, B=( 


/ 





FIG. 2: Changes in the symmetry as the parameters are varied. Our 
starting point denoted by 'SU(4)' corresponds to the choice A — 
B = C^D,E = F = G = Q. 'Duality' V maps a model with 
B > D onto another with B < D and vice versa. Note that the 
full symmetry for generic cases (including the spin-orbital case) is 
SU(2),pi„ X Z2 . 



Of course, (Si — 82)° ± 2i (Si x 82)° (a = x, y, z) give rise 
to singlet-triplet transitions and play a role of raising/lowering 
operators of the pseudo-spins. 
Therefore two building blocks 

^2=^ X/ (Sl,r - S2,r)a (Sl,r+1 - S2,r+l)a 

r a—x,y,z 

'^4 = ^ ^ 4 (Si,r X S2,r)^ (Si,r+1 X S2,r + l)^ 



r a—x^y^z 



behave like J2^r^r+i '^^^ J2^r^r+i ^i'^^ respect to U(l) 
(SO(2)) generated by ^ S,^ except that they have additional 
degeneracy coming from the SU(2)-invariance^''. Of course, 
as far as the pseudo-spin degrees of freedom are concerned, 
the roles of — /i ^ 5^ and J2 ^r^r+i ^re played by the rung 
interaction Ti^ and the rung-rung 4-body term Tiy, respec- 
tively. Although the Hamiltonian B{n2+'H4) + FHe + GHj 
may formally look like a collection of three S ~ 1/2 XXZ 
chains (a = x,y,z) in a finite magnetic field, it is not true 
because — S'f and Si — S2 (a 7^ b) do not commute 
with each other. (In fact, they obey the SU(4) commutation 
relations. See Appendix A.) Nevertheless, an analogy to 
S' = 1/2 XXZ chain is still useful because the model Hamil- 
tonian decomposes into (spin-chirality) pseudo-spin 5 = 1/2 
XXZ part and a (real) spin part. That is, the effective XXZ 
part (H2 + Ti-i + FHe + GHj) has non-zero off-diagonal 
elements only for singlet-triplet transitions (pseudo-spin flip- 
ping), while the magnetic part {Hi and 7^3) gives rise only to 
triplet-triplet transitions. Therefore, we may expect that, if for 
some reasons the spin sector gets gapped and the magnetic dy- 
namics is frozen, the low-energy part of the full dynamics will 
be described by the above effective (pseudospin) XXZ model. 



D. Useful Analogies 

1. Pseudospin- 1/2 model 

Although the Hamiltonian (6) looks complicated, there is 
a useful analogy to a more familiar model-5 = 1/2 XXZ 
model. To understand this, we note that the quantity 

= Si,,-S2,. + J (21) 

formally plays a role of S~ in the 5 = 1/2 XXZ problem. 
Indeed, it is not difficult to show 

[5% (S1-S2)"] =2z(Si XS2)" 

[5% 2(Si XS2)"] =-^(Sl-S2)" (22) 

and 

[5^ (S1+S2)"] = 

[5^Qa6]=0. (23) 

The first two equations imply that Si — S2 and 2 (Si x S2) ro- 
tate like 5^ and 5^, respectively, as is expected from Eq. (11). 



2. spin-1 Base system 

Another formulation of our system facilitates us captur- 
ing the physical meaning of the spin-chirality transformation 
U{0). As a first step, we note that the asymmetric part (Si — S2 
and Si XS2) of the Hamiltonian i?7i2 + -07^4 can be written 
as a hopping term of spin-1 bosons: 

BH2 + DHa 

= {B + D)Y, {bUr+l,a + bUi^,b,,a) 
r.a 

+ {B-D)J2 {blabl+l,a + br+l,abr^a) 

r,a 

(a = X, y, z) , (24) 

where the operator bl ^ creates a spin-1 (triplet) boson with 
spin index a on the site r. Since, in the spin language, the r-th 
rung is occupied either by a singlet or by a triplet, the bosons 
ba should be thought of as hardcore particles. Namely, the 
following local constraints on the boson occupation numbers 
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should be imposed: 



I] =Si,,-S2,. + 3/4 = 5; + l/2 = 0,l . 

a—x,y,z 

(25) 

Note that these particles obey non-standard commutation re- 
lations: 

[br,a, ^^,b] = {Sab{l-ri^) - bl^^br.b} Sr,r' ■ 

One of the greatest merits of this mapping is that the above 
bosons are directly related to the order parameters Si — S2 
(antiferromagnetic) and Si x S2 (p-type nematic) as 



bl 



1 



a,r - 2 (^l,'"~^2,r)°+*(Sl,rXS2^r)° 



x,y,z) . (26) 



The spin-chirality transformation (10) is simply expressed as 
a gauge transformation of these bosons: 



bl 



(27) 



As has been mentioned above. He and H7 have simple in- 
terpretation in terms of an effective spin : 



const 
1 



r 

r ^ 

= G^ri^rir+i - 



^r+1 



const. 



(28) 



(29) 



Last, the symmetric part of the Hamiltonian AHi + CHs can 
be recasted as (see Eq.(5)) 



AHi + CHs 

= ^ [(A + C)T,,-T,+i + 2C (T,-T,+i)' - ACHe 



f const. 
(30) 



In this equation, the spin-1 operators (T^)" = —iSabcbl ^bc,r 
act only on occupied sites. That is, triplet projection opera- 
tors Ptripiet('') on both sides of Tr are implied. From these, 
we can conclude that our ladder Hamiltonian is equivalent to 
an 5 = 1 bosonic t- J-like model on a chain^'^; the devia- 
tion from the self-dual models B — D ^ Q introduces pair 
creation/annihilation processes. A special case of it (C = 
—A < 0), where only the biquadratic interaction exists and 
the system has an enlarged SU(3)-symmetry, was investigated 
by Albertini^^. As has been mentioned in section I, this con- 
nection might be quite useful in realizing the unconventional 
phases in the system of ultra-cold atomic gases-'^. 



III. CONTINUUM LIMIT 



the exact Bethe-ansatz solution'*^, we know that the SU(4) 
point (A ^ B = C = D > 0,E = F = G = 0) is 
gapless and the (conformally invariant) field theory describing 
this massless fixed point was obtained by several authors^*^^**. 
It is given by the level- 1 SU(4) Wess-Zumino-Witten (WZW) 
conformal field theory with central charge c — 3 (for a review 
of WZW- and related models, see for instance Ref. 35). For 
our purpose, it is more convenient to use an equivalent free- 
field description in terms of six real (Majorana) fermions^^. 
The derivation is sketched briefly in Appendix C and the ef- 
fective action corresponding to the fixed point is given by^^ 

3 

Wso(6) = - f E (^R^-^R - ^L^^-^L + XrQ^Xr - X^d^xZ) 

a=l 

(31) 

The above action describes six free massless Majorana 
fermions (^lr,Xlr) i'^ ~ lj2,3) which are equivalent to 
six copies of critical 2D Ising models. In general, SO(6)- 
invariant marginally irrelevant interactions should be added 
to 7iso(6) in order to describe the low-energy physics of the 
SU(4) model (see Appendix C). In what follows, shorthand 

notations ^r l = (Crl' ^rl' ^rl) ^i^^ used. 

Using the quantum equivalence between the level- 1 SU(4) 
WZW model and 7iso(6) (or, the level- 1 SO(6) WZW model), 
we can express all the 15 generators of SU(4) in terms of the 
above six Majorana fermions. In general, we may expect that 
local operators have the following expansions: 



r. lattice 



O + e^'^/^No + e-^^'^/^AA* + {-Ifno ■ (32) 



For the SU(4)-generators, they read 



~ + X{^ + e^-'i^M^ 



(33) 

One can write down X^^^ in terms of six Majorana fermions 
•by using the 6x6 representation of the SU(4) generators (see 
Appendix A for the definition of X^ and Eq. (1)): 

z — ' z 
Sl.L/R = ^xCl/R X ^L/R , S2,L/R = — TTXL/R X XUR ,^ , ^ 

2 2 (34a) 

Gab,LlR = -^^URxl^R (a, 6 =1,2,3). 

The staggered part (n^) is given similarly as 



Si = Cr X : 



S2 = xr X XL 

xio , 



(34b) 



where i? is a regularization-dependent constant. Therefore, 
both the uniform- and the staggered correlation of X"^ are 
written as products of two free fermion propagators and be- 
have like x^^. The second part, whose correlation decays as 
x^"^/^, is more complicated and given by a product of six or- 
der/disorder operators of the underlying Ising models^^. 



A. Field-theory description of SU(4) point 



B. Symmetry operations 



In this section, we develop a low-energy approach to our 
problem starting from the SU(4)-symmetric point. Thanks to 



Aside from the internal SO(6)l x SO(6)r symmetry, various 
symmetry operations keep the fixed-point Hamiltonian (31) 
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invariant. Among them, the followings will play important 
roles: 



by using the following spin-independent SO(2) transforma- 

tion7^L,R(^^L,R)for(e,x")'': 



• Time-reversal: (T) 



L.R 



SR,L ^ ~^L,R ; AR,L * AL,R- 



(35) 



A remark is in order here about T. As is well- 
known, the time-reversal operation is anti-unitary and 
the complex-conjugation must be taken after the trans- 
formation (35) is applied. 

• Translation by 1-site: (Ti.ske) 



SR ^ ~SR > SL ^ ?L 



-Xr 



(36) 



In the field-theory language, Ti-site is nothing but chiral 
symmetry generated by 7^. 



• Interchange of upper- and lower chains: {Vu G Z2) 

xIl ^ Cr,l- (37) 



?R,L ^ Xr,L 



• Site parity: (Ps) 

eR,L(^) ^ eL,R(-2;) 

• Link parity: (Fl) 



xIl{^) ^ Xl,r(-x) (38) 



X^{x) ^ xli-x) 



Xtix) 



Pl 



Pl 



-X^(-x) (39) 



• Ising (or Kramers- Wannier) duality: (si) 



SR,L * SR,L 

Xr 



Xr 



a '^1 a 

Xl — ^ Xl 



(40) 



Cl,r = Cl,r cos ^ - Xl.r sin — 
X^R = ^L,R — + X^R COS — 



(42) 



Strictly speaking, different notations should be used to denote 
the original U{d) defined by Eq. (10) and its continuum ver- 
sion (41,42). However, the meaning is obvious from the con- 
text and we shall use the same notation for the two transfor- 
mations. 

In fact, the fixed -point Hamiltonian (31) is invariant un- 
der an even larger chiral (i.e. left-right independent) version 
SO(2)u X SO(2)l of the above SO(2) and this fact gives us a 
hint to find another set of order parameters. For readers who 
want to know more about the derivation of (41) and (42), we 
give it in the Appendix D. 

The following special cases are worth mentioning: 

• 9 = ^: U{ti/2) = T) reduces to the original 
spin-chirality duality transformation^", where the two 
fermions mix with equal weights. FromEqs. (12a-12g), 
it is clear that models with E ^ Q (no T-odd term Ti^) 
form a closed subset of the full Hamiltonian. 

9 9 = it: This corresponds to the exchange of Si and S2 
(discrete Z2-symmetry 7^12). In the case of Majorana 
fermions, it accompanies a sign change: 



'Xl/r 1 Xl/r ~ Cl/i 



9 ~ 27r: Although this is nothing but an identity 
operation in the original spin- 1/2 -language, Majorana 
fermions change their signs: 



?L/R ~ SL/ 



L/R 



Xl/r 



"Xl/r 



This is already anticipated from what we have for 9 = 
TT. That is, the Z2 -exchange between Si and S2 is not 
realized as a simple Z2 symmetry in terms of Majorana 
fermions^^. Note that this sign inversion does not affect 
the physical operators since they are always written as 
fermion bilinears. 



C. Spin-Chirality Transformation in the Continuum Limit 

In this section, we look for the duality transformation for 
the Majorana fermions. To derive an expression of the 'dual- 
ity' transformation valid in the low-energy limit, it is conve- 
nient to use a generalized version (10): 



U{6)= exp 

rErung 



— + Si^r'S2,r 



In terms of spin operators, U[G) is realized in a non-linear 
manner as we have seen in Eqs. (9a,9b). Interestingly, it is 
realized in the continuum limit as 



U{9)=Tl^{9)Tl^{9) 



D. Order Parameters as Duality Doublets 

In section II, we have pointed out that two quantities 

S1-S2 and 2(SixS2) 

form a doublet under the spin-chirality rotation U(0). How- 
ever, as our system is one-dimensional, the expectation val- 
ues of these vector order parameters are identically zero: 
(Si - S2) = (2(Si X S2)) = 0. Instead, we adopt the fol- 
lowing two rotationally invariant order parameters: 



/-/-) lattice 



(41) 



— Si_r'Si^r+l — S2,r'S2,r+l 

(Sl,r + S2,r) • (Sl,r+1 X S2,r+l) 

+ (Sl,r X S2,r) • (Sl,r+1 + S2,r+l) 



(43a) 
(43b) 
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In numerical studies^', it was shown that phases which are 
characterized by non-vanishing 0^^"'^'^ or Ojjj?"^'^ do exist in 
the phase diagram of the two-leg spin ladder with a ring- 
exchange interaction. The appearance of the latter describes 
an exotic phase since a non-zero value of Osc implies a non- 
magnetic order with T-breaking. 
FromEqs. (11), (43b), and 



it follows that these two order parameters transform as an 
SO(2) doublet: 



/^lattice \ 



cos t/ — sm ( 
sin 6 cos 6 



/^lattice 
/nlattice 



(44) 



That is, as far as the spin-chirality SO(2) is concerned, Osd 
and Osc behave like (Si — S2) and 2(Si x S2), respectively. 
In particular, by the duality V, two order parameters C^p"^" 
and C's^'"^'' are interchanged"'^. 

Now let us find the continuum expressions for the above 
order parameters. As has been mentioned in section III A, any 
local operator on a lattice has an expansion of the following 
form: 

^lattice - Oix)+e'^''^'^Af{x)+e-'^^/^Af*ix) + {-ir0^x) . 

(45) 

By taking operator-product expansions (OPE), we obtain the 

and O^i^'' 



expressions of the staggered parts (C) of O^^^'"^ and ^lattice. 



Osd = * (eR-el - Xr-Xl) and Osc = « (6i-Xl + Xr-^l 

(46) 

By using Eq. (41), it is straightforward to verify 



may recombine them into two scalars and two vectors under 
Ui9): 



scalar: 



vector: 




I Cr-Cl - XR-XL + «(Cr-Xl + XR-6.) 
I (r-Sx - Xr-Xl - ii^R-Xh + Xr-Cl) 



(48) 



(49) 



While the latter two have already appeared, the former are 
new. Therefore, it is suggested that we should add two more 
order parameters, which are scalars under the spin-chirality 
rotation U{9), to complete our analysis. Below, we shall use 
the following set of four order parameters'^: 



CsD = * ( Cr-Cl - Xr-Xl 



Osc = *KR-XL + XR-a 



(50) 



As we already know, the spin-chirality SO(2) transforms 
the first pair (Osd, Osc) as a doublet, while keeping the sec- 
ond (Oq, Orq) invariant: 



0^^i[^r-Cl + Xr-Xl] 
Crq = i Ur-Xl - Xr-Ci 



U{0) 



05 



cos f — sill f 

sin 6 cos 



*-^SD 
*-^SC 



O^ 



(51) 



A similar property holds for the second pair as well. To see 
this, let us introduce the following product'^: 



'-^SD 



smf 



As has been mentioned before, the relations B — D and 
E = define self-dual models, which are invariant under the 
continuous rotation U{6) G SO(2). From the fact that the 
above two order parameters transform as an SO(2)-doublet, it 
readily follows that 



///-jlatticev //-/-)lattice\ ri 

\*-^SD / — \^SC I — ^ 



(47) 



(52) 



which is chiral (i.e. left-right asymmetric) and probably non- 
local in terms of the original lattice spins. Then, it is straight- 
forward to show 



0SD 

05 



f cos 9 
Vsiii^ 



— smt 
cos 6* 



o^ 



(53) 



60 



for generic models on the self-dual manifold' 



E. Second set of order parameters 

It would be interesting to look for the possibility of other or- 
der parameters. Let us restrict ourselves to those which are (i) 
spin singlet (that is, spin indices are contracted), (ii) Lorentz- 
invariant (i.e. Lorentz spin=0) and (iii) have a scaling dimen- 
sion unity. Apparently, we have four such operators ^r-^l, 
Xr'Xl, ^r'XL; and Xr'Cl made up of Majorana bilinears. We 



As in the case of U {9), we will use the notation T) = U{9 ~ 
7r/2) to denote the second duality. The existence of the above 
two dualities (one is non-chiral and the other is chiral) will be 
useful in understanding the low-energy physics and the global 
phase structure. The transformation properties of these order 
parameters under the symmetry operations in section III B are 
summarized in Table I. 

However, this is not the end of the story. For the phases 
denoted by Q and RQ, it will turn out that the q = tt (i.e. 
period-two) components are not sufficient for the full charac- 
terization of the phases. This point will be discussed in section 
IV D in conjunction with the ground-state degeneracy. 
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order param. 






symmetry 






duality 


T 


Tl-silc 


Pl2 


Ps 


Pl 


Sl 


V V 


OsD 


+1 


— 1 


-1 


-1 


+1 




Osc inv. 




-1 


-1 


-1 


-1 


+1 




Oh inv. 




+1 


-1 


+1 


-1 


+1 




inv. OIq 




+1 


-1 


+ 1 


+1 


-1 




inv. Oq 



TABLE I: Four order parameters and discrete symmetries. Note tliat 
under si, D, and D the four order parameters transform onto eaclr 
other. 



interaction as (ge = gj = 0): 

V = .giVi + 32 V2 + .93 V3 + 54 V4 + .95^5 
= ^(51 + 32 + 53 + 94) (Cr-Cl + Xr-Xl) 
+ ^(.91 - .92 + 53 - 54) (6i-Xl - Xr-Cl"! 



+ 2 (-91 + 92 - .93 - 54) - Xr-Xl 



(gi - 52 - 53 + 54) (Cr-Xl + Xr-Cl 



95 ^5, 

(56a) 



F. Interactions in the continuum limit 

Now that we have identified the order parameters, we pro- 
ceed to constructing interactions in the low-energy effective 
action. Basically, we have three different contributions to the 
interactions. The first one comes from the uniform (9 = 0) 
terms in the continuum expressions of the spin operators (see 
Eq. (33)) while the other two from the contraction of the 
±2A;F-terms or of two Akp terms (fcp = 7r/4). 

By using the Majorana expressions (34a,34b) for the spin 
operators, we can rewrite the non-oscillatory part of the seven 
interactions (4a-4g): 





= (Cr-Cl)' + (xr-Xl)' 


+ (Cr-Xl)' + (XR-a)' 


(54a) 


V2" 


= (Cr-Cl)' + (xr-Xl)' 


- (Cr-Xl)' - (XR-a)' 


(54b) 


V3" 


= 2 


(Cr-^l)(xr-Xl) - 


(&-Xl)(xr-^l) 




(54c) 


V4" 


= 2 


(Cr-^l)(xr-Xl) + 


(Cr-Xl)(xr-^l) 




(54d) 


VI 


= 2(^r-|'l - Xr-Xl)(6i-Xl + Xr-Cl) 




(54e) 






^ (^r-Xr + Cl-Xl^ 






(54f) 


V? 




^(6i-Xr)(?l-Xl) • 






(54g) 



The superscripts '0' denote terms coming from the uniform 
(9 = 0) components of the (local) generators. Actually, prod- 
ucts of two n"^ -components contribute the same (but with a 
permuted order) terms to the interactions Vi, . . . , V7, which 
can be absorbed by the redefinition of the bare couphngs. For 
this reason, we will suppress the superscripts '0' and consider 
the interaction 



V = ^5.V. 



(55) 



1=1 



in what follows. These expressions are consistent with the 
transformation properties (12a)-(12g) underW(6') and the dis- 
crete symmetries of the lattice model (see Eqs. (35)-(39)). 
Therefore, we may conclude that with an appropriate choice 
of bare couplings the above seven interactions will describe 
(weak) perturbations from the SU(4)-invariant model. 

It is also very suggestive to rewrite the main part of the 



which can be rewritten into a more convenient form in terms 
of the order parameters; 



V = -^(.91 +.92 +.93 +.94) 

-^(5i- 52 + 53 - 54)(Orq)' 

-^(.91 +.92 -.93 -54) {O^uf 

-^(5i- 52 - 53 + 54)(Os"c)' 



(56b) 



- 55 CsdC'sc • 

From this, we may expect that one of the four competing 
quantum phases is selected according to the values of the 
seven couplings gi , . . . , 177 in the low-energy limit. In order 
to investigate the low-energy behavior of the seven couplings, 
we shall carry out a RG analysis in the next section. 



IV. RG ANALYSIS AND FOUR COMPETING ORDERS 

A. /3-functions 

As has been shown in the last section, we have seven inter- 
actions around the SU(4) (or SO(6)) fixed point. Since we are 
interested in the spontaneous breaking of time-reversal sym- 
metry, we will suppress the TC^ interaction^': = 0. The 
sixth interaction He is a kind of magnetic field or chemical po- 
tential (see section II D) and may be incorporated after physics 
for 96 = is understood. For these reasons, we first consider 
the model with E = F = 0: 



(57) 



At the one-loop level, the calculation of the RG /3-function 
reduces to that of OPE defined below''-: 



Hi{z, z)Hj{w, w) 



\z — w\ 



-nk{w,w) (i, J = 1,2,3,4,7) 
(58) 

The non-zero OPE coefficients listed in Appendix E enable us 
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to write down the RG /3-function^^: 

ffi = 5i + .92 + 5 .93 + .94 
92 = '2 9192 + G 93 94: + 54 57 

53 = 6 51 53 + 2 52 54 (59) 

54 = 2 51 54 + 6 52 53 + 52 57 
57 = -16(5153 - 5254) , 

where dots denote the derivative with respect to the RG time: 
5 = dgi/ (dhiL). 

It is interesting to observe that this set of RG equations 
(RGE) determines a gradient flow in a 5-dimensional space of 
coupling constants. In fact, if we make a change of variables: 

hi= gi , /i2 = .92 , hs = V253 + —^97 , I 



a28 



/l5 = 



128 



57 , 



54 



(60) 



the RGE can be derived as 

from a single RG potential: 

V{hi,h2,h3,h4,h5) = 



dV 
dhi 



1. 



/ll/lo 



hih^ 



hih. 



hihl 



+ Vdhihshs - 3^/2ft.2/^3/l4 - VT0h2hih5 . (61) 

A similar property was pointed out in quasi one-dimensional 
electron systems''^' to explain the simple structure of the phase 
diagrams. 



B. Symmetries of /3-functions 

1. Sign-change- and permutation symmetries 

The above set (59) of /3-functions is invariant under dis- 
crete symmetries (sign change and permutations). To be con- 
crete, the RGE is invariant under the si transformation (see 
Eq. (40)): 

(51, 52,. 93, 54, 0,0, 57) ^ (51, 52, -.93, -54, 0,0, -57) . 

(62) 

Of course, the spin-chirality duality V maps a set of cou- 
plings as 

(51, 52,, 93, 54, 0,0, 57) (51,54,53,52,0,0,57) (63) 

and keeps the self-dual manifold 92 ~ 54 invariant. Combin- 
ing this with si, we obtain for the second duality: 

(51,52,53,54,0,0,57) (51,-54,53,-52,0,0,57) . 

(64) 

For the second duality T), the self-dual manifold is character- 
ized by 



52 



-54 



(65) 



2. Self-dual manifolds 

The existence of the two duality transformations V and T> 
manifests itself in the invariance of the /3-functions under the 
interchange: 52 ^ ±54 (note that the time -reversal breaking 
term is already suppressed). That is, the RG-flow is sym- 
metric with respect to the four-dimensional self-dual mani- 
folds defined by 52 = ±54- Since it follows from the /3- 
function (59) that 

(52 T 54) = (251 T 653 T 57) (52 T 54) , 

systems which are self-dual initially will remain so even af- 
ter renormalization. In fact, the original lattice model with 
52 = 54 has a continuous symmetry U{d)-the spin-chirality 
rotation by an arbitrary angle d-and the Mermin-Wagner- 
Coleman theorem guarantees the condition (72 = 54 is pre- 
served for all orders of perturbation. On this self-dual mani- 
fold, the /3-function simplifies to 



5? 



252 



553' 



51 = 

.92 = .92(251 + 653+57) 
53 = 65153 

57 



25i 



16(51.93 



52 J 



(66) 



This reduced set of /3-functions has a further invariant mani- 
fold characterized by SU(3) symmetry: 

51 = 53 , 

where we have only three coupled equations: 



252 



51 = 651 

52 = 85152 + 5257 
57 = -16(51-52) 



Other invariant manifolds 



(67) 



Furthermore, by adding and subtracting equations in (59), 
we obtain several invariant manifolds: 
Spin-Orbital (SO) manifold: 



91 = 52, 53 =54, 57 = . 



(68) 



This means that as far as we consider the SO model (see 
Eq.(16)), the rung-rung four-body interaction will iiever be 
generated radiatively on manifold. In fact, the underlying 
SU(2)xSU(2) (SO(3)j X SO(3)^) symmetry guarantees that it 
persists even in all orders of perturbation expansion. On this 
manifold, only two couplings gi and 93 suffice for describing 
the low-energy physics and the RG /3-functions substantially 
simplify to the following two coupled equations: 



51 = 25i + 655 



.93 85153 



(69) 



When 93 < and 93 < gi (which may be interpreted as 
the SO-model with J > K/A. For a detailed discussion of the 
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SO-model, see Ref. 57.), the system is attracted in the infrared 
limit by an asymptotic trajectory: 

9*i{^ 9*2) = -.93 -.94) 

(see a region shown by darker gray in Fig. 3) and breaks the 
translational symmetry spontaneously^^. The exact ground 
state'*^ at K = 4 J/3 and the analysis for small K{<^ J)^^ are 
consistent with the above result. The region gi < 0, (73 < \gi\ 
(a portion of Fig. 3 marked by lighter gray) is the basin of the 
gapless SU(4)]^ fixed point. 
Dual Spin-Orbital (dSO) manifold: 

91 = 94, 92 = 93, 97 = (70) 

This is the dual (2?) partner of the above invariant manifold 
obtained by applying V to the SO-model and its existence is 
guaranteed by a hidden (dual) SU(2) x SU(2) symmetry of the 
model. Now two couplings gi and 52 describe the system and, 
again, /^-functions reduce to the set of equations (69) with 53 
being replaced with 32 ■ For 172 < and 52 < <7i, the system 
flows toward another asymptotic trajectory: 

54) = -.92 (= -53) 

and, in this case, breaks time-reversal symmetry T (see the 
next subsection for the detail). The basin of SU(4)]^ fixed point 
of the SO manifold is now mapped onto gi < 0, 172 < In 
terms of lattice models, this corresponds to the model Hamil- 
tonian 

with J < K. Similarly, we have two more invariant manifolds 

91 = -52 , .93 = -.94 , .97 = , ^^^^ 
91 = -94 , 92 = -93 , ,97 = , 

on which the full RGB reduces to two coupled equations. Four 
RG-invariant lines (rays) defined by the intersections among 
these invariant manifolds (hyperplanes) will turn out to corre- 
spond to four dominant competing phases. 

Away from these invariant manifolds, the rung-rung four- 
body interaction gj will be generated in the course of renor- 
malization even though it is absent initially. 



^3 SU(4) FP 




FIG. 3: RG flow for the SO- and dSO model. The origin (shown 
as SU(4) FP) corresponds to the level- 1 SU(4) WZW model and a 
line gi — 33 to self-dual models (SU(4)-symmetric, in this case). A 
dashed line roughly corresponds to a path across the SU(4) Suther- 
land model (extended massless— >SU(4)^staggered dimer (SD)) dis- 
cussed in literature (e.g. Ref. 57). An asymptotic trajectory in 
the lower-right portion (highlighted by a thick arrow) is responsi- 
ble for the spontaneous breakdown of translational- (for SO model) 
and time-reversal (for dSO model) symmetries. 



where the constant to marks the crossover point where the 
weak-coupling perturbation breaks down. Plugging these into 
the RG equations and requiring consistency, we have a set 
of non-hnear equations which determines the IR-asymptotics. 
We found various solutions and among them the followings 
are relevant for our analysis: 
Self- dual: 

(Gi,G2,G3,G4,G7)= Q,0,-i 0,^^ (I) (73a) 
(Gi,G2,G3,G4,G7)= (^'O>^'O,-0 (11). 

(73b) 

The meaning of these rays will become clear if we plug Eqs. 
(73a,73b) into (56b): 



C. Asymptotic form of RG-flow 

The set of RG /^-functions (59) is complicated and com- 
plete analysis of it is not so easy. We numerically integrated 
the equations and found that the RG flow exhibited a striking 
feature in the low-energy limit. To see this more closely, we 
apply an ansatz proposed by Lin, Balents, and Fisher'^ in the 
study of a half-filled two-leg Hubbard ladder 

Since only marginal interactions appear in the /3-functions, 
a natural ansatz valid for the infrared asymptotics may be'^ 

5.W = T^ (z = 1,2,3,4,7), (72) 
to — t 



(1): - ig; [{Ol^f + {Olcf) + p^Hr (74a) 
(II): - ^ffil (to' + W) - ^.9n^7 . (74b) 

Therefore, the above two rays (1) and (11) may be thought of as 
corresponding to SD-SC- and Q-RQ transitions, respectively. 

From the argument in section 11, it is obvious that the latter 
case-(ll) has SU(3) symmetry, while the SU(3) symmetry of 
the former is hidden and appears only after the particle-hole 
transformation for the left movers is applied (see section V A). 
Also it is worth noting that they are related to each other by 
the Ising duality si. 
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S0( 6)-symmetric rays: 



D. Four dominant phases 



(Gi, G2, G3, G4, Gj) 



(Gi, G2, G3, G4, G7) 



(Gi, G2, G3, G4, G7) 



(Gi, G2, G3, G4, G7) 



11 1 1 

8' 8'"8'"8' 



1 1 



1111 



1 1 1 



(SD) 

(75a) 

(SC) 
(75b) 

(Q) 

(75c) 

(RQ). 
(75d) 



Note that the latter two are invariant both under spin-chirality 
SO(2) and under SU(3), while the former two are not. It is 
important to note that the low-energy effective Hamiltonians 
for these rays assume the following simple forms: 

Ha = Hso(6) - g*A (Olf (A = Q, RQ, SD, and SC) . 

(76) 

Semiclassical argument based on Eq. (56b) tells us that 
these symmetric rays correspond to the four competing phases 
where one of the order parameters has a non-zero expectation 
value: {0\) ^ 0. In particular, from Table 1, it follows that a 
spin-disordered ground state with broken T is realized along 
the ray 'SC . Furthermore, we can show that all these rays cor- 
respond to IR-massive flows along which translational sym- 
metry is broken and the SO(6) symmetry is restored asymp- 
totically (see next section). 

Transitions among dominant phases: 



1 



Gi - G2 = - , G3 = G4 = , G7 = (Q 



SD) 



(77a) 



Gi 



-G2 



1 



Gi — G4 



Gi 



-Ga 



, G3 = G4 = , G7 = (RQ ^ SC) 

(77b) 

, G2 = G3 = , G7 = (Q <^ SC) 

(77c) 

' G3 = , G7 = (RQ ^ SD) 

(77d) 



None of these restores the original SO(6) symmetry. Actu- 
ally, these rays correspond to transitions among the above four 
dominant phases. 

Case of four competing orders: 



Gi — 1 , G2 — G3 — G4 



G7 







(78) 



This ray describes non-trivial quantum criticality resulting 
from the competition among four different orders, which will 
be discussed in section VI. 

The point here is that in the low-energy limit our system 
is asymptotically characterized only by a single (diverging) 
coupling constant'^'^-' ~ l/(^o — t)- 



1. SO(6)-restoration 

To illustrate how SO(6) symmetry, which is explicitly bro- 
ken in the bare action (57), is restored along the four rays (i.e. 
'SD', 'SC, 'Q', and 'RQ'), we take the two rays 'Q' and 'SD'. 
From the continuum expressions (56a) and (56b), it is obvious 
that the model along the ray 'Q' is explicitly SO(6) invariant: 



Ti-Q ~ Hsotb) — g* (Cq) 
= Wso(6) + 5*(?R'Cl 



Xr-Xl)^ 



(79) 



This is nothing but the Hamiltonian of the SO(6) Gross-Neveu 
model^"*. This model is integrable. Its spectrum is known^^'^^ 
to consist of the fundamental fermion with mass AI together 
with a kink of mass m = AI/^/2. The Z2-symmetry Oq ^ 
— Oq is broken spontaneously and, as a result, we have a finite 
expectation value (Oq) 7^ 0. 

The application of the Ising duality si (Eq. (40)) to this 
Hamiltonian Hq takes us to another ray 'SD' and the asso- 
ciated Hamiltonian Hsd since 7iso(6) Wso(6) and Oq i-^ 
OsD- A similar argument applies to the other pair (TYrq and 
Use) as wefl. 

Moreover, the spin-chirality duality V (V) maps Hsd (Wq) 
onto Hsc CHrq). These facts imply that all these four low- 
energy Hamiltonians Hq, Hrq, Hsd, and Hsc can be trans- 
formed back to the SO(6)-invariant one (79) by applying sym- 
metry operations of 0(6)l x 0(6)r. This is the reason why we 
named the asymptotic rays (75a-75d) 'SO(6)-symmetric'. The 
relationship among these four phases is summarized in Fig. 4. 



Gaussian 




Ist-ordgr 



Ist-order 



Kramers-Wlannier ^^tTtT^ Kramel-s-Wannier 

(si)i (l:l:l:l>SU(4)^l:-l:l:-l)^si) 




Gaussian 
(dual) 



FIG. 4: Relationship among the four dominant phases. Ratio of the 
couplings gi, 32, gs, and g4 along each symmetric ray is shown. 
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2. Ground-state degeneracy 

Now we know that four competing (gapped) phases exist 
around the SU(4)-point and that they are mapped onto each 
other by the two discrete transformations V and V. From this, 
we may conclude that all these phases have the same ground- 
state degeneracy (from the exact solution""' for the SD phase, 
it might seem natural that we should have the degeneracy 2 
for the others as well.) But there is a pitfall here; the dis- 
crete transformation T) that relates 'SD' to 'Q' is not a local 
one and there is no reason to believe that T) maps a phase 
onto another one with the same ground-state degeneracy. A 
typical example of this kind of non-local transformation is 
the so-called Kennedy-Tasaki transformation^^, which maps 
a unique (infinite-volume) ground state onto four-fold degen- 
erate ones. 

To find the correct answer, we first note that Dirac fermions 
used to construct our SU(4) Hamiltonian (see Appendix C) 
are expressed as exponentials of boson operators and that 
the semiclassical bosonic ground states should be considered 
modulo the following gauge-equivalence: 

4>a,a,VR ~ 0a,(TX/R + V^^ Na,c,,LIR {Na.aX/R G , (80) 

which keeps the bosonic expressions of the fermions (B2) un- 
changed. 

We carefully counted the number of inequivalent ground 
states for all four cases 'Q', 'SD', 'SC, and 'RQ'. Although 
the calculation was slightly more complicated than that in Ref. 
15, the procedure was essentially the same (see Appendix F 
for details) and the results are summarized as (see also Table 
11): 

'Q'-phase: 4 ground states; 2fcF-part of (X^X^j^) is non- 
vanishing (a kind of quadrumerized phase) and period-4 
ground states appear 

'SD' -phase: 2 ground states; 2fcF-part is vanishing (staggered 

phase) and the ground states have a period-2 

'SC'-phase: 2 ground states; 2fcF-part is vanishing 

'RQ' -phase: 4 ground states; 2fcF-part is non- vanishing 

(period-4 density wave of the pseudo-spin (Si,r-S2,r))- 

For some reasons, the period-4 {q ~ 7r/2) contribution in 
the order-parameter correlations vanishes for the SD- and SC 
phases, as is expected from the exact ground state^^ for which 
we have {OrOr+i) — (— 1)''3/4 (when the amplitude of the 
order parameter is maximal). On the other hand, we have the 
2fcF-components in the Q- and RQ phases. This implies that 
the 4fcF-component is not sufficient to fully characterize the 
latter two phases. 

So far, the nature of the ground states appearing in the 'Q'- 
or the 'RQ' phase is not so clear. This is mainly because 
the transformation si relating to Oq is non-local and 
we do not know the microscopic (or, lattice) expressions for 
Oqj^q. By symmetry arguments, we can restrict the possible 



forms of lattice operators corresponding to Oq . We have five 
building blocks for the lattice order parameters all of which 
have (at least for the leading contributions) the same 4A:F(or 
q = 7r)-term Oq in the continuum. Exactly on the SU(4)- 
invariant line, these five will condense with the same ampli- 
tude to describe a phase with SU(4) quadrumerization (note 
that we have 4-fold degeneracy). Unfortunately, our contin- 
uum approach cannot tell anything about which of the five 
will become dominant when we move away from the SU(4) 
manifold. 

Similarly, since the 4fcF-part of Si_r-S2,r is written as O^q, 
the RQ-phase may be thought of as a density-wave phase of 
the pseudo-spin 5^. By using the analogy in section 11 D 2, 
this is nothing but the bosonic charge-density wave state with 
period-4. 



phase 




G.S. degeneracy and symmetry 




degen. 


T 


Vl2 


Ps 


Pl 


SD 


2 


+ 1 


-1 


-1 


+ 1 


SC 


2 


-1 


-1 


-1 


+ 1 


Q 


4 


+ 1 


+ 1 


-1 


+ 1 


RQ 


4 


+ 1 


+ 1 


+1 


-1 



TABLE II: Four dominant phases and the ground-state degeneracy. 
In all four phases, translation symmetry is broken. 



V. PHYSICS OF SELF-DUAL MODELS AND PHASE 
TRANSITIONS 

Now that we have identified four dominant phases in our 
SO(6) problem, the next question to ask would be the quantum 
phase transitions among them. In this section, we discuss the 
nature of the phase transition between SD and SC phases as 
well as the Q-RQ transition. As described in section IV C, 
the transition should belong to the self-dual manifolds (172 ~ 
±34) and we first present a fully quantum description of the 
low-energy physics along these manifolds. 



A. SU(3) X U(l) bosonization 

According to the analysis presented in section IV C, the 
one-loop RG flow on the self-dual manifold, which is invari- 
ant under the spin-chirality rotation (6*) or U (9), is attracted 
to one of the following two special rays: (1): gi = —173 = 
91 /G, 92 = 34 = 0, .97 = 4.91* /9 and (11): gi = .93 = .gfj/G, 
52 = 54 = 0, 97 = -4gi*i/9 with g^^^ > 0. Along these rays, 
the low-energy field theory reads as follows: 
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9L 
6 



-1/(11) n, . 9ll 

rtiR = nso(6) + Y 



■4l - Xr-Xl 



Cr-Cl + Xr-Xl 



Cr-Xl + xr-Cl 



Cr-Xl - Xr-Cl 



^ (^r-Xr) (Cl-Xl) 



(81a) 
(81b) 



It is interesting to observe that the two low-energy field the- 
ories Tif^ and Tlf^^ interchange themselves under the Ising 
duality symmetry si (Eq.(40)). 

As has been shown in section IV C, the low-energy effective 
Hamiltonians (8 1 a,8 lb) describe the competition between two 
different orders (see Eqs.(73a), (73b)): the staggered dimer- 
ized and the staggered scalar chirality O'^q orders for the 
model Tif^ (81a), and the two quadrumerized phases Oq and 

CJq for the model TiJ"* (81b). The main question concerns 
the nature of the quantum phase transitions that result from 
the competition between these orders (Ogo ^ O^q for the 
case 1 and Oq ^ O^q for the case 11). In this section, we 
present a full quantum description to this issue. 

The key observation, which will be crucial in the follow- 
ing analysis, is that the low-energy theories (81a,81b) display 
a hidden symmetry which is not SU(2)xU(l) but, in fact, a 
larger U(3) symmetry. In the following, we are going to dis- 
cuss only the case of the first model (81a), i.e. the nature of 
the phase transition between the staggered dimerized- and the 
scalar chirality phases. The physical properties of the second 
model (81b) can then be derived readily by applying the Ising 
duality symmetry si (40). 

In the following analysis, it would be convenient to com- 
bine three pairs of Majorana fields (^^/^ and Xl/r to form three 
Dirac fermions ^'ax/R- 

v&„,R = ^i±M , vI.,,L = ^^^, (82) 



V2 



V2 ' 



with a ~ 1,2,3. The reason for this left-right asymmetric 
definition is that SU(3)-symmetry is clear in this notation. The 
symmetric notation will be useful in describing the second ray 
(II) (Eq. (73b) or the model (81b))''^ 

From Eq. (82), it is straightforward to express the order 
parameters (50) in terms of these Dirac fermions: 



OgQ and O^Q look like the density operators. If we had 
adopted the left-right symmetric definition instead of Eq. (82), 
Ogjj and Ogp would have behaved like the 'Cooper pairs'. 
This analogy can also be revealed by mentioning the expres- 
sion of the spin-chirality rotation Z-/ (6*) (Eq. (41)) and the Ising 
duality symmetry si (Eq. (40)) on these Dirac fermions: 



*a,R 




*a,R , 






^-a.R 


u{e) ,- s 
> e 2 


*a,R , 




*a,R 


^a,R 







u(e) 



^-ax (84a) 

e-*^ *a,L (84b) 
(84c) 



The appearance of the left-right asymmetric expressions for 
U{9) is just an artifact of our definition (82); if the symmet- 
ric definition had been used, Eqs. (84a) and (84b) would 
have been interchanged. One thus observes that the spin- 
chirality rotation acts as a pseudo-charge U(l) symmetry on 
these fermions. 

The next step is to express the low-energy field theory cor- 
responding to model (81a) in such a way that U(3)-symmetry 
of our problem is manifest. To this end, it is useful to in- 
troduce the chiral SU(3)i currents built from the three Dirac 
fermions: 



Jr\=^^' 



a.b 



(85) 



where the 3x3-matrices T^, A — 1, are the generators of 
SU(3) in the fundamental representation 3, which are normal- 
ized so that Tr(r'^r^) = 5^^/2. Using the identity: 



X! '^ab - ^ (SadSbc - ■::SabScdj , 



(86) 



and the definition (82), one obtains the following identity: 



a=l 
3 

a=l 
3 

^^SD = *E(*--R<L + <R*aX 

3 

^^SC = E (*-.R*I,L - *iR*aX 



a=l 



(83) 



In this notation, two order parameters Oq and Orq, which 
are invariant under U (9), transform like 'Cooper pairs', while 



1 



{^R<L 



Xr-Xl 
1 

^ 6 



I + (Cr-Xl + Xr-^l. 
(Cr-Xr) (Cl-Xl) . (87) 



We thus deduce that the low-energy field theory (81a) associ- 
ated to the ray (I) exhibits an exact U(3)-symmetry: 



a=l 



^E^^^^^^- (88) 
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Noting the well-known fact^^ that the three massless Dirac 
fermions (the first term) can be bosonized in terms of a sin- 
gle scalar field and the above SU(3)i currents, we can further 
recast the above Hamiltonian as 



''•IR — 



A=l 



A=l 



(89) 



where tp and d are the U(l) free bosonic field and its dual, re- 
spectively. There is a 'spin-charge separation' between the 
'charge' sector described by the Tomonaga-Luttinger (TL) 
model (first term in Eq. (89)) and the SU(3) non-Abelian 
'spin' sector described by the second- and third term-the 
SU(3) Gross-Neveu (GN) model. The appearance of SU(3) 
is not so surprising; as has been mentioned in section II, we 
can embed SU(3)-symmetry into the triplet (spin-1) sector of 
the original ladder models. Therefore, we may conclude that 
the SU(3)-sector describes the dynamics of the (real) spin de- 
grees of freedom. The physical meaning of the remaining 
'charge' sector will be clarified below by using a slightly dif- 
ferent bosonization scheme. 

As is known from the exact solution^^, for g* > a spec- 
tral gap (spin gap) is formed in the 'spin' (SU(3)) sector and 
the low-energy physics is dictated by massive SU(3) spinons 
and antispinons which transform like the fundamental repre- 
sentations 3 and 3, respectively. The low-energy field theory 
(89) displays nevertheless a c = 1 quantum criticality due 
to the decoupled 'charge' degrees of freedom which are de- 
scribed by a massless free boson field (TL model) at the free- 
fermion point. Therefore, the quantum phase transition be- 
tween the staggered dimerized phase and the scalar chirality 
phase generically continuous and belongs to the c = 1 univer- 
sality class. 

To clarify the role of the 'charge' U(l) ((p) and the 
'spin' SU(3)-sector, we apply the Abelian bosonization to the 
Hamiltonian (88). As in Appendix B, three bosonic fields are 
introduced to bosonize the three Dirac fermions as follows: 



exp (^iv47r^a,R 

CXp ( -iVlTT^aX 



(90) 



Since we have adopted the left-right asymmetric definition in 
the above equations to bosonize our three Dirac fermions, the 
spin-chirality rotation simply reads (see Eq.(84a)) 



Po.L/R 



'Ja,L/R 



(a = 1,2,3) 



(91) 



The next step is to switch to a new basis where the U(l)- 
and the SU(3) degrees of freedom are separated from each 



other: 



ip = 

ipf : 



1 
1 



62 + ' 



+ 02 - 2(/>3) 



(92) 



The two bosonic fields ips and ipf are compactified fields with 
special radii i?s_f so as to capture the underlying SU(3) sym- 
metry of the problem: 

ip,,t - (^s,f + 27ri?s,f, (93) 

with R, = l/^/2^ and Rf = ^3/27r. The spin-chirality 
transformation affects only the 'charge' sector: 



(94a) 




while the dual spin-chirality Z^/ (6*) changes only the dual field 

i9: 



ip ip 




(94b) 



Now the physical meaning of the 'charge' fields (p and 
is clear; the field ip (respectively i9) describes the angular 
(or, phase) fluctuation of the doublet (Og^, Ogc) (respectively 
(Oq, Orq)). In the bosonization approach, two fields ip and -d 
are conjugate (or, dual) to each other Therefore, two doublets 
(C'fD,e'^c) and {0^,01q) are mutuafly dual objects which 
are analogous to the Cooper pairs and the density-wave oper- 
ators in strongly-correlated electron systems. 

Using this Abelian bosonization, it is then possible to write 
down the SU(3) current-current interaction in Eq. (88) in 
terms of the bosonic fields: 



A 



2ti 



cos 



27r 



ips^ COS ^VGtt ipf^ 
+ dx^PiR dxipih) 



(95) 



cos 



A straightforward semiclassical analysis of the potential 
part of Eq. (95) together with the identification (93) show 
that the ground state in the SU(3) "spin" sector is three-fold 
degenerate with expectation values: 



(96) 
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Now the role of the current-current interaction ^ J^J^ is 
clear. For the semiclassical vacuum configurations (iy9s) and 
((/7f), the modulus of the order-parameter doublet (Csdj '^sc) 
is given by 



6 



X^A ^R^L- 'Th^ oiily remaining massless fluctuations in the 
azimuthal direction are described by the TL (or, Gaussian) 
model (the first term in Eq. (89)). 



That is, fluctuations of the doublet in the radial direc- 
tion is suppressed by the marginaUy relevant interaction 



J 



It is then interesting to express the order parameters of Eq. 
(83) in terms of these bosonic fields; 



1 

TTflo 
1 

TTflo 



2 COS [yi 
2 cos ( V2Tr 



1 



IT i9s) COS (V47r/3 + i?f) + cos (^^/Air/i i9 - ^/Sir/i i?f) 

Vs^ COS 7r/3 i^f^ + cos (^\/ 47r/3 iy9 — ^Stt/S (^9/ 

2 cos ^\/27r lySs^ sin ^■\/47r/3 ip + \/2t:/3 ipfj + sin (^y/ 47r/3 — ^^871/3 (^f^ 
i?s^ sin 7r/3 i9f 



2 cos I V 27r 



(97a) 
(97b) 
(97c) 
(97d) 



r 



It is straightforward to observe that for the type-(I) self-dual 
model, i.e. for the model described by the field theory (81a), 
one has (Og) = (OJd) = (OJc) = (O^^) = due to the 
quantum-criticality of the "charge" degrees field in Eq. (89). 
From Eqs. (97a-97d), we also deduce that the first doublet 

and O^Q has correlation functions decaying as 
i.e. has quasi long-range coherence, whereas the second one 
Oq and O^Q is exponentially decaying due to strong quan- 
tum fluctuations. The situation is completely reversed for the 
second (type-(II)) self-dual model. 



B. Small deviation from selfdual models 

Now let us discuss the effect of small deviation from the 
self-dual model (81a). The deviation may be incorporated by 
adding the following symmetry-breaking perturbation to the 
Hamiltonian (89): 



^SB — £ (^2 — 'H4) 



^r-Xl + Xr-Cl 



(98) 



We can also express this symmetry-breaking perturbation in 
terms of the three Dirac fermions using Eq. (82): 

Use = 2e^ [*a,R*l,L*b,R*L + "^Ir^ a.L^l^^bX 

a.b 

(99) 

The effect of this term can be elucidated by means of the 
Abelian bosonization (90) of the Dirac fermions. Moving to 
the basis (92), we obtain the bosonized form of the symmetry- 



breaking term (98) 



27r (^s) cos (2y^4Tr/3 ip - ^2^ ipf'^ 



+ cos (2y^ATr/3 ip + ^/Stt/S ipf'j 



(100) 



The stability of the critical line described by the model (89) 
with respect to the small perturbation (100) with |e| ^ 1 can 
be investigated by a naive semiclassical analysis: for |e| ^ 1, 
the bosonic fields ip^j are still frozen to one of the ground-state 
configurations ( (<^s) , (ipf) ) and the non-abelian "spin" degrees 
of freedom are still gapfull. 

The only difference from the previous case is that here 
we have couplings between (p and (/3s,f and they may shift 
the (semiclassical) values {{^Ps), iff)) from those obtained for 
B = D (Eq. (96)). However, as far as the value of e is small 
enough, we may expect that (/jj.f still are pinned at the same 
values. Assuming that (^s_f are locked to, for instance, the first 
set in (96), we deduce that the U(l) "charge" degrees of free- 
dom acquire now an extra term: 



4e 



Hc^- [d^'pf + {d^-df -^^cos 2V4V3¥' 
2 L J tt^Oq V 

(101) 

The low-energy field theory for the "charge" degrees of free- 
dom takes thus the form of a quantum sine-Gordon model. 
The interaction has scaling dimension A = 4/3 < 2 so that 
the perturbation is relevant and the "charge" degrees of free- 
dom acquire a gap. For e > i.e. (72 > 54 (e < i.e. 
92 < 94), the U(l) bosonic field is locked on: ((^) = or 
{ip) = \/^/2 ({ip) = \/3^/4 or {ip) = 3^3^/4). Using the 
identifications (97a-97d), we then deduce: (Osd) 7^ 0, and 
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(©5) = (^^sc) = (O^o) = for e > 0, i.e. for 32 > 54, 



SO that one enters the staggered dimerized phase SD. In con- 
trast, when e < 0, i.e. 32 < 94, we have: (O^q) 7^ 0, 
and (05) = (OJd) = (Crq) = and the scalar chiral- 
ity (SC) phase is stabilized by the small symmetry-breaking 
term. A similar result can be obtained by considering the sec- 
ond ground state (ips) = ^7r/2, (i^f) = In Fig. 5, we 
illustrate how the fluctuation of the order parameter doublet is 
frozen as the symmetry is lowered. 

Finally, we remark that a similar approach can be applicable 
also to the competition between the quadrumerized phase (Q) 
and the RQ phase described by the low-energy field theory 
(81b). Repeating the same steps as before, we see that the 
theory obtained by making the replacement tp {) in Eq. 
(101) dictates the competition. 





(a) 



(b) 




FIG. 5: Schematic picture illustrating how the doublet of order pa- 
rameters is fixed: (a) At the SU(4) point, the doublet fluctuate both 
in the radial- and angular directions, (b) For the self-dual models, 
the modulus of it is fixed by the interactions in the SU(3)-sector, 
whereas fluctuations in the azimuthal direction are still gapless. (c) 
The remaining angular fluctuation gets locked by a small deviation 
from the self-dual models. 



C. On the possibility of umklapp 

In the previous subsection, we have seen that on the self- 
dual manifold the enlarged U(l)-symmetry forbids the cosine 
interactions with explicit (^-dependence from appearing in the 
low-energy Hamiltonian. However, we are on a lattice and 
any interactions compatible with both the U( 1 ) symmetry and 
the discrete lattice symmetries are allowed in principle. Now 
let us discuss briefly about the possibility of a gap generation 
by such umklapp interactions. 

To this end, it is necessary to know how the discrete 1-site 
translation is realized in terms of our bosons. Apparently, the 
1-site translation concerns only the charge boson ip and its 
dual i9. In order for the expression of the translation (see Eq. 
(36)) 

(Cl, XL, £,R, Xr) {£,L, XL, -61, -Xr) 

written in terms of the Majorana sextet to be translated cor- 
rectly to a bosonic expression, we should take 



(102) 



From Eqs. (97a-97d), it is easy to check that it indeed changes 
the sign of four order parameters 0\ (note our order parame- 
ters are staggered ones). Using this, we can obtain a selection 
rule for the umklapp. 

In what follows, we shall look for terms which (i) does not 
contain (/3-field and are (ii) translationally invariant. To this 
end, we note that the charge bosons (p and 1) enter into the 
expressions of (chiral) Dirac fermions 'I'a.L/R like 



4'a,R 



where the ellipses denote the contributions from (p^ (t?s) and 
ipf (df). In particular, 'Cooper pairs' are U(l)-invariant and 
will take the following form: 



Since an operator 



(*a,R*6.L 



acquires a phase Nn after the 1-site translation, the integer N 
should be even. This gives us the first constraint. 

However, this is not the whole story. In general, interac- 
tions constructed this way will contain not only (p^ and ipf but 
also their duals i?s,f which are not pinned by the interaction 

^r3^l'^ if these interactions include the dual fields they 
will be suppressed by strong quantum fluctuations (as was 
pointed out by Schulz^*' in the context of spin chains). For this 
reason, we have to look for combinations with (i) N =even 
and (ii) no explicit i?s k&i) dependence (these interactions will 
be generated by perturbations even if they do not exist in the 
bare Hamiltonian). By a direct enumeration, we checked no 
such interactions appear up to = 4. For = 6, we have 
several combinations e.g. 



(«'i,R*1,l) (^'2,R^'2,l) (^'2,R^'i,l) 
(4'2,R*3,l) (4'3,R^2,l) (^3,R^3,l) 



(103) 



(note that due to the Fermi statistics these expressions should 
be understood as short-distance expansions). In general, terms 
containing both ^'a^R^'tx and ^ are also allowed. 

However, up to A^ = 6 we did not find such combinations 
satisfying our requirements. (In principle, we can consider in- 
teractions made up of A^/2 pieces of 'I'a.R^'h.L and the same 
number of ^\ ^. But they do not contain ?9-field and are 
not umklapp.) 

Therefore, possible umklapp interactions with the lowest 
dimensions will be of the following form: 



cos 6 




const. 



(104) 



Around the SU(4) point, it has scaling dimensions larger than 
two and is thus irrelevant; strong renormalization of the Lut- 
tinger K is needed for the umklapp (104) to stabilize a gapped 
phase on the self-dual manifold. 
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VI. OTHER QUANTUM CRITICALITIES 

Now let us discuss the properties of the low-energy effec- 
tive Hamiltonian on the asymptotic rays (77a)-(77d) and (78) 
which correspond to the transitions among the four domi- 
nant phases. Combination of the Abelian- and non-Abelian 
bosonization enables us to obtain non-perturbative solutions 
to the problem. Although it looks more complicated to treat 
the last one (78), technically it is slightly simpler and we will 
begin with the case of four competing orders (78). 



A. Case of four competing orders 



On the ray (78), the effective action reads as follows: 



in the SU(3)-part, we may naively expect additional critical 
degrees of freedom resulting from this symmetry. However, 
this is not the case. To see this, we first note that the interac- 
tion part of Eq. (105) is written compactly as 



Wint = -y (^Cr X ?R + XR X XR ) • ( X + XL X XL 



) 

(107) 



where we have introduced the level-2 SO(3) currents Ir^l de- 
fined by; 



^ (6l,L X 61, 



L + Xr,l X xr,l , 



(108) 



(Osc) 



y E (^R^-^R - ^L^.e^ + Xr^.Xr - Xld.xl) 



Hq-SD-SC-RQ = ^50(6) - Y (^(Csd) 

+ (Oq)' + (Orq)^ 



Cr • Cl + (xr • Xl) 



^R • XL 



Xr • ^1 



(105) 



with g* > 0. 

As has been mentioned in section IV C, the model (105) 
describes the competition between the four different orders of 
our problem: the staggered dimerized (SD), scalar chiral (SC), 
and the two quadrumerized (i.e. period-4) orders (Q and RQ). 

We apply the SU(3)xU(l) bosonization scheme to our 
Hamiltonian (105). After some algebra, we obtain the fol- 
lowing simple Hamiltonian; 



"^Q-SD-SC-RQ 



^ E ["^"^ a; • ) + COS ( VStt di 



i=l 



(106) 



where = {if^,ipf), 6 ~ (??s,?9f), and ai are the posi- 
tive roots of the SU(3) algebra: di = (1/2,^/3/2), as = 
(l/2,-^/3/2), and as = (1,0). Eq. (106) describes a 
Lie-algebraic generalization^"'^' of the self-dual sine-Gordon 
models considered in Ref. 72. 

From Eq. (106), one immediately observes that Hamilto- 
nian (105) displays U(l) quantum critical behavior due to the 
non-interacting bosonic field ip. The interacting part of (105) 
takes the form of an SU(3) self-dual sine-Gordon model with a 
marginal interaction. From the self-duality symmetry (f. 



It is interesting to note that Xl,r are simply written as 

Ir,l = ^U(5")^R,L , (109) 

where {S°')hc — ~i£abc (a = x, y, z) is a 3x3 representation 
of spin operators which form an SO(3) subgroup of SU(3). 

By using these SO(3)-currents, we can then rewrite the ini- 
tial Hamiltonian (105), which governs the competition of the 
four orders, in a current-current form: 



Wq-sd-sc-rq =^ (^{dx^f + (dx'&f^ 



(110) 



The second part of this Hamiltonian is that of the level-2 
SO(3) WZW model perturbed by an SO(3)-invariant current- 
current interaction and is integrable^"*. Now the relation to 
the previous model (Eq. (89)) is clear; if one adds current 
interactions coming from the remaining five SU(3)-currents, 
the 6'-dependent part of (106) is canceled and the form (95) is 
recovered. 

For g* > 0, which is relevant to our problem, the second 
part has a spectral gap and a non-trivial structure of massive 
spinon. We thus deduce that the initial Hamiltonian (105) ex- 
hibits only a U(l) quantum criticality due to the free boson (p. 
The transition that results from the competition between the 
four orders is thus of a U(l) Gaussian type. 



B. Transitions among dominant phases 

The transitions between two dominant phases within the 
same (spin-chirality) doublet (i.e. SD-SC and Q-RQ) have 
been akeady discussed in the previous section in conjunction 
with the low-energy physics of the selfdual models. Here we 
concentrate on the transitions between different doublets (e.g. 
Q-SD). 

Let us consider, for instance, the Q-SD transition (Eq. 77a) 
and begin with rewriting the interaction term (Oq)^ + (Osd)^- 
9 Plugging the bosonized expressions (97a) and (97b) into the 
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above, we obtain 
1 



1 ^ 

+ ^ [cos(V87r/3r-'^) + cos(V87r/3r-i?) 

TT Oq L 

(111) 

where ~ {Lp,(ps,(pf), i? = i?s, i?/), and the last 
summation here is taken over all six positive roots of 
SU(4): = (l/2,±V3/2,0), (1/2, ±1/(2V3), T^aTs), 
(0, 1/V3, v/273), and (1,0,0) normalized so that \Pr\ = 1- 
Therefore, the effective Hamiltonian corresponding to the ray 
(77a) reads 



T-Cq sd - IT 



* 6 



|^cos(%/87r/l-(^) + cos(%/87r/3r-??) 



(112) 

where we have rescaled the velocity v so that {dx0Y terms 
coming from interactions may be absorbed. If we replace 
SU(3) (a) by SU(4) (/3) and (/) by in the previous section, 
this is nothing but the second part of (106). Therefore, the 
argument here goes similarly to that in the previous section; 
the effective Hamiltonian (112) can be written here in terms 
of SO(4)-currents J^]^ (1 < i < j < 4) as''^ 



Q-SD 



i]\2 



(113) 



2g' 



Kj 



That is, we have obtained the level-2 SO(4) WZW model 
(with central charge c = 3) perturbed by an SO(4)-invariant 
marginal current-current interaction which is an integrable 
field theory^"*. In contrast to the previous case, we have no 
critical degrees of freedom here and we expect that the transi- 
tion between Q-phase and SD-phase is of first order A similar 
result holds for RQ-SC transition as well. 



VII. GLOBAL STRUCTURE OF THE PHASE DIAGRAM 

In this section, we shall investigate the main effects of the 
interactions (Hsfi) in Eq. (6) that we have so far neglected 
in our field-theoretical approach. In addition, we shall also 
use variational and strong-coupling analyses to figure out the 
global phase diagram of model (6). 

A. Effects of 

So far, we have neglected the T-breaking three-spin inter- 
action 7^5 of Eq. (4e). In this section, we shall give a hand- 



waving argument for its main effect. In section V, it has been 
argued that the T-breaking and the appearance of the SD/SC 
phases can be understood as a pinning of the vector doublet 
(OJp, OJc) the spin-chirality plane (see Fig. 5). Now the 
role of H5 can be discussed in a similar manner by means of a 
semiclassical argument applied to the bosonized Hamiltonian. 

According to Eq. (56b), the T-breaking interaction 
changes the form of the effective potential as 



9_ 
6 



^SD 

6 



(O^sd) + {Old 



(114) 



SsC'sdCsc 



In the presence of Tis, the spin-chirality plane is no longer 
isotropic and the principal axes of the ellipsoid (see Fig. 5) 
are tilted by 7r/4 (the sign of 55 determines the direction of 
the longer axis.). Then, we may expect that the two order pa- 
rameters and O^Q simultaneously take non-zero values. 

To make the above argument more quantitative, we add the 
gs-term to the interaction (95). The interaction takes its min- 
ima when ((/3s) and ((/3f) are given by Eq. (96). Then the ex- 
pectation values of the cp-field are determined by finding the 
minima of trigonometric functions. For {tp^) = {(pf) = 0, (ip) 
is determined by minimizing —9.95/2 sin(-\/l67r/3(y5). The 
calculation goes similarly for the other cases as well. 

This analysis tells us that the doublet (Osdi ^sc) takes one 
of the following values when T-breaking interaction 55 is 
present; 



(0SD,C)SC) 




for 55 > 
for 95 < 



(115) 



The above argument can be easily generalized to the case of 
non-self-dual models. 



B. Effects of Tie 

Now let us discuss the effect of the He interaction (Eq. (4f)) 
which has been neglected in the preceding analysis. As has 
been mentioned in section 11, our model (spin) Hamiltonian 
can be rewritten in terms of spin-1 (hardcore) boson &r,a (a = 
x,y, z) as: 

n = AUi + BH2 + CHa + DHi + EH5 + FH& + GH7 

= (i? + i^) ^ (6,t A+l,a + A,a' 



-E'^Eabc (bl^abl+i^bKc + &J,c^r,a^^r+l,6 + (r ^ r + 1' 

r 

+ ^ [(A + C)T,-T,+i + 2C {Tr-Tr+if] 



G^n>^+i+ -4C + F--G^<, (116) 



21 



where the projection onto occupied states is implied for the 
fourth term (the triplet-triplet interaction). From this, it can 
be easily read off that the coupling F of He affects the chemi- 
cal potential of the hardcore boson. Therefore, for sufficiently 
large values of the system becomes either a carrierless in- 
sulator (i.e. ^ Q for all sites r) or a fully occupied state 
(i.e. — 1 for all r). Apparently, the ground state of the car- 
rierless insulator is trivial and is, in terms of the original spins, 
given by a tensor product of local (rung) singlets. The fate of 
the spin wavefunction in the latter case depends strongly on 
the two couplings A and C which dictate the spin-spin inter- 
action between hardcore bosons. Here, we will mainly focus 
on the case A + C > 2\C\ which corresponds to the Haldane 
phase^''. 

When |F| is decreased, we have quantum phase transi- 
tions to a conducting state. Let us first consider a transition 
from the carrierless insulator Then, the transition is equiv- 
alent to the superfluid(SF)-onset transition of spin-1 bosons. 
For the moment, let us neglect the interactions among dif- 
ferent species of particles (including the hardcore repulsion). 
As is well-known^^, the bosonic two-body interactions are 
relevant at the z = 2 SF-onset transition and the system 
flows toward the strong-coupling fixed point. In the pres- 
ence of U(l)-breaking anisotropy (B ^ D), the particle num- 
ber is no longer conserved and the fixed point is replaced 
by that of the z = 1 critical 2D Ising moder^-^**. Since 
we have three copies of such Ising models, we may expect 
that the fixed point here is described by the level-2 SU(2) 
WZW model which is equivalent to three massless Majorana 
fermions with c — 3/2. Now let us estimate the impact of the 
neglected boson-boson interactions. Since they are given by 
0(3)-invariant (i.e. spin-symmetric) products of four Fermi 
fields, the only possible one should be the marginal interac- 
tion of the form Jr Jl (Jl,r are the level-2 SU(2)-currents 
of the WZW model). Therefore, we may conclude that the 
SF-onset transition from the rung-singlet phase {carrierless 
insulator) to the staggered dimer- (SD) or the scalar-chirality 
(SC) phase belongs to the level-2 SU(2) WZW criticality; the 
SD- or SC order appears as a consequence of the SF ordering 
in the presence of the U(l)-breaking anisotropy. This is con- 
sistent with the results of weak-coupling approach^^'^^'^*^ and 
numerical analysis^^. Since the number of the bosonic parti- 
cles is conserved, the transition at B = D (self-dual) is ex- 
ceptional and belongs essentially to the Japaridze-Nersesyan- 
Pokrovsky-Talapov (JNPT) universality class^^ (note that we 
have three copies of the JNPT criticality here). 

A similar argument applies to the transition from the fully- 
occupied state as well. The elementary particle here is a 
rung singlet moving in the sea of rung triplets (note that we 
have assumed that we are in the spin-gapped phase and the 
magnon excitations occur at the energy scale of the 'Haldane' 
gap). Then the dynamics of these excited singlets may be 
modeled by the interacting spinless fermion. The sign of the 
fermion-fermion interaction will be positive (negative) when 
G is much larger (smaller) than A and C. When A and C 
are not so large and the anisotropy is absent {B = D), the 
F-driven transition is of the JNPT universality class again. 
The xy-anisotropy (B ^ D) alters the transition to that of the 



Ising-type. We show a schematic phase diagram illustrating 
the effects of Ti.^ and Tig in Fig. 6. 
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FIG. 6: A schematic picture illustrating the effects of Tie and the 
anisotropy e — B — Din the spin-chirality plane on the global struc- 
ture of the phase diagram of our model. According to the sign of the 
anisotropy, we have two different orderings; an anisotropic SF in the 
x-direction (T-even) for e > and one in the j/-direction (T-odd) 
for e < 0. The duality symmetry D interchanges the region with 
e > and that with e < 0. On both sides of the region enclosed by 
a dashed line (dominant singlet fluctuations), we find conventional 
phases (spin-1 or rung-singlet), which can be interpreted as insulat- 
ing ones. 



C. XXZ-analogy and pseudo-spin description 

In order to demonstrate the relevance of the pseudo-spin 
picture (see Section II D 1), we present two approaches which 
are complementary to that developed in the previous sections. 
In both approaches, the key is how to identify the pseudo-spin 
degrees of freedom on the lattice. 



1. Variational approach 

Let us begin with the self-dual models. In the above, we 
have seen how the U(l) Gaussian model emerges as the low- 
energy effective field theory after the spin sector is gapped. 
The spin and the (pseudo)charge decouple from each other at 
low energies and the physics of the spin sector may be de- 
scribed by the following bilinear-biquadratic interaction (see 
section II and Appendix A 2 for more details) 



AHi + CHs 



E 



[{A + C)T,-T,+i + 2C(T,-T,+i)2] 



For this reason, we may expect that the Haldane-gap physics a 
la Affleck-Kennedy-Lieb-Tasaki^'' dominates in a reasonably 
large region of the phase diagram. Of course, we have the 
additional pseudo-charge degrees of freedom (i.e. motion of 
triplet rungs in the rung-singlet background) here and the sta- 
bility of the Haldane state is not so obvious. However, from a 
simple argument'*'^ we know that the valence-bond-solid state 
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proposed in Ref. 75 is stable against the motion of singlet 
rungs (or, holes in the spin-1 t-J model given by Eqs. (24)- 
(30) or Eq. (116)). In fact, numerical simulations carried 
out for the spin-1 bosonic t-J model**' (the set of parameters 
B = D = t/2, A = J, C = Q was used there) show the 
existence of a finite spin gap and the Luttinger-liquid behav- 
ior in the pseudo-charge sector, which is consistent with our 
conclusion in section V A. 

Now that we know that the spin sector is generically gapped 
and decoupled from the pseudo-charge sector, the next ques- 
tion would be what kind of degrees of freedom determines 
the global structure of the phase diagram. From the bosoniza- 
tion analysis presented above, a natural guess would be the 
pseudo-spin 5=1/2 (or 'charge' in the t-J language). As 
was described in section II D, 7^2,4 correspond to the pseudo- 
spin-flipping processes and Tig and Tiy to magnetic field and 
the 5^5^-interaction, respectively. The remaining parts con- 
cern dynamics in the (true)spin sector, which is separated from 
the low-energy sector by a finite spin gap. 

To develop a variational theory for this kind of spin liquid, 
we first construct a coherent state of the pseudo-spin S by 
combining a local singlet |s),- and a triplet \i)r: 



rGrung 



COS 



2 



tan ■ 



It),. 



where 0(r) and 6{r) are respectively azimuthal- and polar an- 
gles of the spin vector at rung r. Apparently, this is unsat- 
isfactory because we still have local spin degrees of freedom 
represented by |t) r. To kill the spin degrees of freedom and 
construct a spin-gapped wave function out of the above co- 
herent states, the most natural way would be to replace triplet 
states \i)r on r-th rung by the following 2x2-matrix'^^''^'': 



|t0>r 

V3 \+^\t- 



l/r 



-V2\ti) 

-|t0>r. 



(\ta) is the rung-triplet state with T^ = a) and use the modified 
ansatz: 



<S> f\/3e*'^('')tan^|s),.l 



rGi'ung 



(117) 



Note that this is essentially the same as the building block of 
the matrix-product ground state adopted in Ref. 46. Taking 
the trace of the above matrix product, we can obtain the de- 
sired (unnormalized) spin-singlet wave function with a finite 
correlation length. This wave function may be thought of as 
the valence-bond-solid (VBS) state^^ randomly 'diluted' by 
vacancies (see Fig. 7(b)). 

To illustrate the usefulness of our pseudo-spin picture, we 
map out the variational phase diagram of the two-leg spin 
ladder with four-spin cyclic exchange 7Yiadder+4-spin (see Eq. 
(20)). As a variational ansatz, we choose 



{4>i,9i) for r=even 
((/)2, O2) for r=odd . 



(118) 



Thanks to the special form of jfi), we can reduce the computa- 
tion of the ground-state expectation values (f^|7Yiaddei+4-spin|f^) 
to that of 4x4 matrices^^-^-' and we obtain the variational en- 
ergy i?var(0i, Oi,(t>2, O2) to minimize. Since the explicit form 
of Ey3i-{(f>i,9i,4>2,02) is unimportant, we only show the re- 
sulting phase diagram in Fig. 7. It should be compared with 




1 Vs ^4 



self-dual point 



(b) 



FIG. 7: (a) Variational phase diagram obtained by the ansatz (117). 
Shown is the 'pseudo-spin' direction (tf), 9) for each phase: con- 
ventional rung-singlet- and rung-triplet (Haldane or VBS) phases 
are spin-polarized states while the SD- and SC phases may be 
viewed as canted states with different XY-projections. (b) Typical 
spin-singlet configuration contained in our variational wave function 
(117). Ovals and open circles respectively denote triplet- ('occupied 
sites') and singlet ('unoccupied sites') rungs. The spin-1 bosons to- 
tally form a spin-singlet valence-bond state. Due to motion and pair 
creation/annihilation of triplet bosons, the wave function (117) de- 
scribes a strongly fluctuating state. 

the phase diagram of the same model obtained by large-scale 
numerical simulations^^ (DMRG). The path searched in Ref. 
29 corresponds to the line Jr = 1 in Fig. 7. Three dominant 
phases (rung-singlet, SD, and SC) found there appear in our 
variational phase diagram as well, while the phase denoted 
by 'dominant vector chirality'^^'-^'^, whose nature is not clear 
currently, is missing in ours. 

Similarly, the ferromagnetic phase is beyond the scope of 
our simple variational calculation since our variational wave 
function (117) contains only spin-singlet phases. In our pic- 
ture, both rung-triplet- and rung-singlet phases correspond to 
(pseudo)spiyi-polarized states as is shown in Fig. 7. The spin- 
chirality transformation U{9) rotates the pseudo-spin vector 
along the z-axis. When the external magnetic field (Tie) is not 
very strong, spin-canted states (SD and SC) are favored by the 
xy-coupling Ti2,4; in the SD (SC) phase, the xy-projection of 
the pseudo-spin lies in the x(y)-direction. 



2. Strong-coupling approach 

Another way to see the role of the pseudo-spin degrees of 
freedom would be a strong-coupling expansion. A natural 
starting point might be isolated rungs (as in the usual lad- 
der systems). Unfortunately, however, all the phases that we 
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have found in the field-theory analysis break the periodicity 
of the original Hamiltonian and the limit of isolated rungs is 
not quite helpful (the phases accessible from the limit would 
be more or less trivial). Instead, we divide the whole lattice 
into plaquettes (see Fig. 8(a)) and introduce a new coupling A 
which controls the coupling between neighboring plaquettes. 
On these plaquettes, a kind of coarse-grained degrees of free- 
dom will be defined by which the low-energy sector of our 
ladder system may be described. 

To explore the vicinity of the SU(4) point, we adopt the 
following Hamiltonian: 

J2 Jr{hi{r) + il + 52)h.2{r) + il + d3)h3ir) 

r— rung 

+ (1 + Si)hiir)} +5eY, Si,r-S2,r , 



(119) 



where the couplings Jr alternate like 



Jr = I (for r even) , J,, = 1) (for ?- odd) 

(120) 

and hi (r) denote local Hamiltonians obtained by restricting 
Hi to two rungs r and r + 1. 
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FIG. 8: (a): A two-leg ladder as a weakly-coupled plaquettes. 
Strongly-coupled plaquettes, on which pseudo-spin l/2s and coarse- 
grained spin- Is are defined, are shown by thick lines. All interac- 
tions (both two-spin and four-spin) connecting these plaquettes are 
multiplied by A(^ 1). (b) and (c): Two triplet states which con- 
stitute a (nearly) degenerate ground-state manifold in the vicinity of 
the SU(4) point. Two spin-l/2s contained in each oval form a spin 
triplet. In both (b) and (c), two triplets (ovals) are anti-symmetrized 
to form a total spin triplet on a plaquette. 

Let us begin with the case of A = 0, i.e. the limit of isolated 
plaquettes. On a single plaquette, we have 16 states which can 
be decomposed into the following multiple ts: 



(0)^0(1)^0(2) 



(121) 



Around the SU(4) point, two of the three triplets (see Fig. 
8(b,c)) form nearly degenerate ground states (exactly at the 
SU(4) point, they are completely degenerate). The explicit 
forms of these states are given as 



^ (122a) 



l[]>^;^(l|| ID-Ill ll»^;^(0 + |->)- 

(122b) 

where arrows and double lines in the above respectively de- 
note spin singlets and triplets (ovals in Fig. 8). 

It should be remarked that we shall construct an effective 
theory in terms of new degrees of freedom (spin-1 and species 
of triplets) defined on these plaquettes. In this plaquette pic- 
ture, our ladder is already dimerized and the two staggered 
phases (SD and SC) should be understood as uniform ones. 
Correspondingly, we shall introduce order parameters defined 
on each plaquette: 



Osd{p) — Si,r'Si^r+l ~ S2,r'S2,r+l 

Osc(p) = (Si,. X S2,,) • (Si,,,,+i + S2,.+i) + (r 



M), 
(123) 



where p labels the plaquette formed by four points (l,r), 
(2,r), (2,r + l),and (l,r + 1) (see Fig. 8). 

The point here is that the eigenstates of Osd and Osc can 
be constructed out of the two triplets shown in Fig. 8: 



I ± 1)SD 



V2 



(124) 



where ±1 denote the eigenvalues of C'sd,sc(p)- 

The next step is to identify the states having the pseudo- 
spin (r) up and down. To this end, we consider the coherent 
state of a spin- 1/2: 



in) = cos -It) 



e+*"^ sin ■ 



i) 



If we identify the following eigenstates of and 

k^-±i/2) = ^(IT)±U)) 

|r^-±l/2) = ^(|T)±^U)) 



(125) 



with the eigenstates ( 1 24) of the order parameters (Osd and 
Osc), it is suggested that we should take 



T)-|[]) 



i) 



!) 



and 



Osd(p)=2t; , Osc(p) = 2rp^ 

Tp = Si,r-S2,r + Si,r+1-S2,r+1 ■ 



(126a) 



(126b) 



Obviously, is a generalization of discussed in section 
II D to a plaquette system. 

Now that we have identified the pseudo-spin degrees of 
freedom, it is straightforward to carry out a perturbation ex- 
pansion for (nearly) degenerate ground-state manifold made 
up of 1^^) and J| ). Within this manifold, operators for 
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spin, p-nematic, and n-nematic read 

Sl,r + S2,r = + 1 + S2,r+1 = ^ 1 ® Tj, (127a) 

Si^r ^ S2,r = ^ S2,r+1 — T'^ <E) Tp (127b) 

2(Si,, X S2,.) = 2(Si^,+i X S2,.+i) = ® Tp (127c) 
n"'3 — ( qa qP I qP qo' \ 



the pseudo-spins: 



(127d) 

where Tp denotes a spin-1 operator on the plaquette p. Plug- 
ging these equations into the original Hamiltonian, we ob- 
tain the (Ist-order) effective Hamiltonian which describes the 
physics around the SU(4) point: 

Heff = 

A ^ |(1 + <52)r;r;+i + (1 + S^y^^rl^, + U Tp-Tp+i 

p— plaq ^ 

+ A(l + <53) {t;t;^,) {Tp-Tp+i + 2 (Tp-Tp+i)' - 2} 

p 

+ ^ ((52 - (53 + (54 + (56) Tp^ 

p 

= A ^ |(l + 52)r;r;+i + (l + (54)rp%V + J}tp-Tp+i 

p=plaq 

+ 2A(1 + (53) (r;r;^,) \ ^ QfQ^+i - 1 



(128) 



where we have introduced the coarse-grained n-type ne- 
matic operators Q^'^ = (TpT^ + T^Tp)/2 on plaquettes. 
This is similar to the well-known Kugel-Khomskii effective 
Hamiltonian"*^ for the orbital-degenerate systems; here the 
role of two degenerate orbitals is played by two types of 
triplets I ) and | [ J| ) ■ It should be noted that the xy- 
components of the pseudo spin (AF and p-nematic r^) 
couple to the dipolar part (Tp Tp+i) of the magnetic Hamil- 
tonian while the z-component is associated with the (n-type) 
nematic part (Q^^'Q"^). 

In order to further simplify the effective Hamiltonian, we 
again assume that the spin sector is well described by the VBS 
wave function^^. Then the dipolar- and the nematic part of 
the effective Hamiltonian may be replaced by the following 
expectation values: 

(Tp.Tp+i)vBS--| , (Y.QfQp+i) 

/vBS 

(129) 

and the above effective Hamiltonian (128) reduces to that for 



^A Y {(l + '52)r^rpV + (l + (54)rp%V + j} 

p— plaq 



-A(l + (53) ^Tp%V 

p 

+ Y(^2-S3 + S4 + Se) Tp" 



(130) 



This is nothing but the S* = 1/2 XYZ Hamiltonian in a mag- 
netic field^''. For small 5^ the system is xy-like (easy plane) 
and ferromagnetic ordering will occur mainly in the xy-plane; 
when 62 > S4 {62 < 54) the pseudo spins align in the x{y)- 
direction and the SD (SC) phase will be stabilized (see Eq. 
(126b)). For 62 ~ S4 (the self-dual models), on the other 
hand, the system reduces to the well-known XXZ Hamilto- 
nian and the low-energy sector is generically described by the 
TL model (B6b). These are exactly what we have found by a 
field-theoretical analysis in section V. 

For appropriate choices of S2, S3, and 84, the zz-part will 
dominate and staggered ordering along the z-axis will occur 
(note that the sign of r^r" is antiferromagnetic). Since this 
ground state has a period two in the 'plaquette' picture, it has, 
when translated back to the original ladder, a period 4. 



VIII. CONCLUDING REMARKS 

In the present paper, we have developed a unifying ap- 
proach to the problem of unconventional phases appearing in 
generalized two-leg spin ladders. Although the low-energy 
properties of the two-leg spin ladder with only two-spin (or, 
ordinary exchange) interactions are fairly well-known, not 
much is known for the case with four-spin interactions. 

To treat the problem from a wider viewpoint, we have con- 
structed, out of S* = 1/2 spins and their bilinears, a general 
lattice model Hamiltonian with SU(2) (rotational) invariance 
and symmetry under the exchange of two constituent chains. 
We have adopted the so-called SU(4)-symmetric model as the 
starting point to elucidate the interplay between two-spin- and 
four-spin interactions. In particular, the SU(4) symmetry is 
the maximal continuous symmetry that one can have in gen- 
eral two-leg 5=1/2 spin ladders. The greatest advantage of 
this extended symmetry approach is that it can unify the un- 
conventional phases (e.g. SC phase) of the two-leg spin ladder 
with 4-spin exchange interactions. 

In the present paper, we have shown that four unconven- 
tional phases (Q, SD, SC, and RQ) of the generalized two-leg 
spin ladder are unified at the SU(4) multicritical point. These 
phases, which emerge in regions characterized by substantial 
strength of the four-spin interactions, is quite difficult to de- 
scribe by means of a perturbative approach starting from the 
limit of two decoupled chains. Fortunately, however, we have 
a well-controlled starting point-the SU(4) point-in our prob- 
lem. In the field-theory language, the low-energy properties 
of the SU(4) point are described by the level- 1 SO(6) WZW 
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model which is equivalent to six copies of the 2D Ising model 
at its critical point. The continuum expressions of the inter- 
actions can then be written in terms of this basis and we have 
revealed the existence of a duality symmetry V, the so-called 
spin-chirality duality^", and a new hidden one T) which is an 
emergent symmetry. The one-loop RG flow exhibits quite a 
simple structure (four dominant phases and symmetric rays 
toward them) in spite of the complexity in the RG equations. 
The duaUty transformation V and its dual (V) map the dom- 
inant phases onto each other (see Fig. 4). The existence of 
these duality symmetries enables us to investigate the nature 
of the quantum phase transitions between these competing or- 
ders which occur along the self-dual manifolds. The behavior 
on the self-dual manifold B = D is intriguing; the spin sector 
dies away by opening a gap and in turn the pseudo-spin de- 
grees of freedom come into play in the low-energy physics. 
Their gapless spin-singlet fluctuations are described by the 
c = 1 TL model. In fact, what controls the competition and 
quantum phase transitions between the unconventional spin- 
singlet phases (SD and SC) is this pseudo-spin degrees of free- 
dom and we demonstrated this picture by a simple variational 
calculation in section VII C. The fact that the pseudo-spin 
sector displays an U(l) quantum critical behavior helps us un- 
derstand the global phase structure in the same manner as we 
do in the S* = 1 /2 XXZ model. Note that the above U( 1) crit- 
icality is not an emergent one, whereas the other one within 
the manifold 92 = —.94 is emergent in the sense of Ref. 14. 

Another interesting point concerning the physics on the 
self-dual manifold is the emergent SU(3) symmetry. Although 
this is not so obvious in the lattice analysis, our field-theory 
analysis teUs us that instead of the original SO(3) an en- 
larged SU(3) symmetry appears in the low-energy limit. As 
a byproduct, we established an interesting description of our 
problem in terms of the spin-1 (hardcore) bosons. In this re- 
spect, we hope that our results obtained for spin systems will 
be of some help in understanding the phases of the spin 1 bo- 
son systems . 
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APPENDIX A: SUMMARY OF SU(4) 
1. Generators 

In section II, we introduced 15 SU(4) generators in the 
form of 4x4 matrices. Below, we give the explicit forms of 
them. Note that we use the singlet-triplet basis defined by 
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(A2b) 

The remaining 9 generators are essentially given by tensor 
product of two Pauli matrices; 



ia,b^ 1,2,3) . 
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(A2d) 



(A2e) 
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2. SU(3)xU(l) 

The spin-chirality transformation introduces a natural cat- 
egorization of the above generators in terms of a subgroup 
SU(3)xU(l). The easiest way of understanding the appear- 



ance of SU(3) would be to decompose the whole Hilbert space 
into a singlet and a triplet (of course, the singlet corresponds to 
U(l)-factor); the SU(3)-part acts only to the triplet subspace. 

With this identification, the SU(3) Gell-Mann matrices are 
given as follows: 
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They are supplemented by the generator of (spin-chiraUty) 
U(l): 



G, 



U(l) 



(Gil + G22 + G33) 



which satisfies 

[Gu(i),G'']=G (a = l,...,8). 

These 9 operators generate a subgroup SU(3)xU(l). 

It is a straightforward task to rewrite our building blocks 
Eqs. (4a-4g) in terms of the above generators. In particular. 
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one has: 



AHi + CH^ 



0Q = —(j)a,R + 0a. L- We impose the following commutation 
relation for the chiral bosonic fields 



= AiG^{r)G^{r+l)+G''{r)G''{r+l)+G\r)G'^{r+l) 
+ c{G\r)G^{r + 1) + G^{r)G^{r + 1) + G\r)G\r + 1) 
+G6(r)G6(r+l)+G«(r)G*^(r+l)+Gu(i)(r)Gu(i)(r+l) 



(B4) 



so that the left- and the right movers of the same species anti- 
commute: {^'a.R(x), 4'q,l(2/)} = 0- The anticommutation 
between fermions of different species is guaranteed by the 
presence of the Klein factors (here Majorana fermions) in 
the definition (B2), which obey the anticommutation rule: 



(A3) 



{ka, iib} — 2Sab- 



(B5) 



From which we observe that the model with A = C has an 
SU(3) symmetry. 

Interestingly, the remaining six SU(4) generators 

S1-S2 , 2(81x82) 

are order parameters of the above SU(3)xU(l). In fact, 

(81 -82) + i 2 (81x82) and (81 - 82) - i 2 (81 x 82) 

respectively transform like 3 and 3 under SU(3). The spin- 
chirality U( 1 ) acts like 

(81 - 82) ± i 2 (81 X 82) ^ c^'" [(81 - 82) ± i 2 (81 X 82)] . 

APPENDIX B: BOSONIZATION DICTIONARY 

To estabhsh the notations, we briefly summarize the main 
formulas used in our bosonization analysis. Consider free 
Dirac fermions defined by the following Hamiltonian: 

a=l 

(Bl) 

Let us first introduce N bosonic fields to bosonize N Dirac 
fermions as follows"*^: 



exp 



(a = l,...,iV) 



exp -iv47r0QX 



(B2) 



flo being some ultraviolet cut-off (typically the lattice spacing) 
and we use the following normalization for the chiral bosonic 
fields: 



K,R (z) 06,R (w)) = - ^ In (z - w) 
47r 

K,h (z) (t>b.h (w)) = --2- In (2 - w) 
All 



(B3) 



where the complex coordinates z and z are defined by z = 
f r+i X and z — vr—ix, respectively. The total bosonic fields 
4>a and their dual fields 6a are defined as: = (f'a r + 4>a l. 



Using these bosons, TYoirac can be written as a Hamiltonian of 
the so-called Tomonaga-Luttinger (TL) modeP^: 



(B6a) 



fdxY,{Jla + Jla) 

•' a=l 
a=l 

(B6b) 

where the operators of the form should always be un- 
derstood as normal-ordered '. A" '. The above equiva- 
lence is established by the following bosonization formulas 
for fermion bilinears: 
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(B7a) 



*l,R*a,L = 
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An 

(B7b) 

J^I^R^x^-a.R (B7C) 

(B7d) 



/47r<I>„ 



APPENDIX C: SU(4) EFFECTIVE ACTION FROM 
HUBBARD MODEL 

In this section, we sketch the derivation of the SU(4) fixed- 
point Hamiltonian from the two-band Hubbard model at quar- 
ter filling: 

"^Hubbard = ^ ^ (^'^li,a,i+l ^a,a,i + C^a,n,i Ca,(T,j+l^ 



U 



Y^ ria,a,inb,a',i{^ " Sa,bSa,a') , (CI) 



i,a,b. 
(7, a' 



where ^ a creates an electron with spin cr(=t, [) in orbital 
a(= 1,2) of the site-i and Ua.a.i = c\^^ iCa,a,i- The key 
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idea is the following. When the Hubbard interaction U is 
sufficiently large, the system becomes a Mott insulator and 
the charge degrees of freedom gets decoupled from the low- 
energy physics^^'^^. However, we still have four states (a = 
1,2, cr =t, I) on each site and we identify them with four 
states of the fundamental representation 4 of SU(4). These 
states constitute the low-energy sector of the system. As in 
the usual single-band case, the effective Hamiltonian describ- 
ing the low-energy physics is obtained by the second-order 
perturbation: 



a, a 

b, a' 



4^2 

IT 



E 



Pi.i+l 



4*2 
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^XfXr+i+ const. , 



(C2) 



where the projection onto the subspace with one electron per 
site is enforced by Pg. The operator Pi.i+i appearing in the 
last line is the SU(4) permutation operator acting on 4(8)4 on 
neighboring sites. The Hamiltonian ^ Xf is nothing 
but the SU(4)-invariant model (17) in section II. 

If we repeat the same steps starting from the continuum ex- 
pression of the two-band Hubbard model (CI), we may ob- 
tain the desired continuum field theory for the SU(4) model 

(17)5<-58. 

Let us begin with the first term of Tieubbard Eq. (CI). As 
is well-known in the standard bosonization treatment of one- 
dimensional electron systems, the spectrum of the four elec- 
trons (ci = ci_|, C2 = ci j, C3 = C2,-f, C4 = C24) near the 
Fermi points k = ±kp = ±7r/4 may be linearized to ob- 
tain the Hamiltonian of the Dirac fermions (B 1 ) with iV = 4 
and V ~ 2t sin kp. In order to obtain the continuum expres- 
sion of the second term (Hubbard interaction Hu), we need 
the expression of the density operator rib i. Since we have 
two species of Dirac fermions (L and R) for each fermion, the 
electron density is made up of four pieces; 



ao 



(C3) 



(6=1,2,3,4) 



With the help of Eqs. (B7a)-(B7d), the right-hand side can be 
readily bosonized as: 



Using this bosonized expression, one finds that the effec- 
tive Hamiltonian which describes the low-energy properties 
of the two-band Hubbard model at quarter-filling consists of 
four terms: 
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The two parts 'H2kr and 7i4fep contains terms oscillating with 
wave vectors q ~ ±2kp and q = ikp respectively and may 
be neglected in the first approximation. Before presenting the 
explicit forms of He and 7Ysu(4), it is convenient to move to 
another set of boson basis: 



(C6) 



Since the above is an orthogonal transformation, the Dirac 
part is transformed to 



Hoirac = ^ ^ { (9, $a)' + (c'^Oa)' } . (C7) 

a— c,s,f,sf 



The [/-dependent part (Tic + 'Hsis(4)) can be recasted as: 
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U 



cos(v47r$s) cos(v47r$f) 
cos(V^$f) cos(V47r<&sf) 
+ cos(V47r$sf) cos(\/47r$s) 



where M is a 4x4 matrix with all matrix elements equal to 
unity. From Eqs. (C7) and (C8), we can read off 
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(C4) 



dx \ ^ + (9.ec)') + ^id-^A (C9a) 
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and 

'WsU(4) 



'' a=s,f,sf ^ 

H — ^ /rfx jcos(V47r$s) cos(%/47r$f) 
+ cos(A/47r$f) cos(\/47r$sf) 



cos(v47r$sf) cos(v47r$s) 



(C9b) 

which is written only in terms of three bosons $s, $f, and 
$sf and dictates the physics of the SU(4)-sector Comparing 
with the bosonization treatment of the single-band Hubbard 
model, one finds that the (Akp) umklapp scattering present 
in the single-band case is absent here while the SU(4)-part 
('SU(4)' should be replaced with 'spin' in the single-band 
case) contains the marginally irrelevant backscattering as in 
the case of the single band. This may seem contradicting since 
we know the two-band Hubbard model will become insulat- 
ing (i.e. charge-gapped) for sufficiently large U {U ^ 2.8t)^^. 
This paradox is remedied by taking into account the second- 
order (U^) perturbation coming from the OPE Ti.ikf'HAkf In 
fact, the above OPE yields the contribution like 

3[/2 



which will pin the charge $c-field for large enough U . There- 
fore, we may fix the value of the 4>c-field and drop Ti^ in the 
following analysis. 

The final step is to rewrite the SU(4) Hamiltonian 7isu(4) in 
terms of six Majorana fermions. The mapping is provided by 
the following correspondence; 



XivL + i Xr/l 
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cxp(-ipV47r(^s,R/L) 
cxp(+zpV47ri^t;R/L) 
cxp(-ip%/47r.^st\R/L) , 



(CIO) 



where the integer p is defined by p = +1(L) and p — — 1(R). 
The Majorana fermions ?/a (a = 1,2,3) have been intro- 
duced to guarantee the anticommutation. The following for- 
mulas proved by taking OPEs are useful in rewriting the SU(4) 
Hamiltonian (C9b): 



dx^lRll. = «%/7rXR/LXl/L 
dx^sfm. = *%/^^R/LXRrt. 

cos(A/47r$s) = i7rao(Ki - 
cos(V47r$f) = iTTao{K4 - 
cos(\/47r $sf) = z7rao(/v3 
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The Majorana fermions are normalized so that 



2n{z - w) (C13) 

(rwx'H) = 0, 

and the Majorana bilinears are defined as 

Ka=^^^^ , Ka+3=XRXt {a ^1,2,3). (C14) 

By using Eqs. (Cll) and (C12), the first term in Eq. (C9b) 
can be written as: 

3 „ 



(C15) 



U dx (kiK2 + + KsKq) , 



where we have introduced the renormalized velocity v — {v — 
U /{2tt)). Similarly, Eq. (C12) enables us to recast the second 
term of 7^50(4) as: 



-■^ j dx {^^2/^<}j +U Jdx{KlK2 + 



K4K5 + KsKe) . 

(C16) 

From these equations, we finally obtain the desired result: 

- 9^ E / (^R^-^R - ^L^-^L + Xr^.Xr - xldxxD 



(C17) 



For repulsive interaction U > 0, the interaction (the last 
term) is marginally irrelevant and known to yield logarithmic 
corrections^^'^^ to physical quantities. 



APPENDIX D: DERIVATION OF THE DUALITY 
TRANSFORMATION IN THE CONTINUUM LIMIT 

If the duality transformation (10) works not only on a lat- 
tice but also in a low-energy field theory, the uniform part of 
the generators derived above should obey the same transfor- 
mation rule (since the duality transformation does not change 
momentum of the system). Plugging Eqs. (42),(34a)into (11), 
we obtain the following set of equations for the uniform part 
ofSi: 

rl" = 

cos' (0 ee + sm' (0 X^x'-^sin^rx'-eV) , 

(Dl) 
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where we have dropped the indices L/R and {a,b) = 
(2, 3), (3, 1), (1, 2). Since all three equations are satisfied si- 
multaneously for a transformation which is independent both 
of color (a = 1, 2, 3) and of the chirality (L/R), we try the 
following: 



(D2) 



After some algebra, we see that the above equations are satis- 
fied if we choose the parameter as 




APPENDIX E: NON-ZERO OPE COEFFICIENTS 



The operator-product expansion of interaction operators is 
defined as: 



'Hi{z, z)'Hj{w, w) 



\z — w\ 



.Hkiw.w) (i J = 1,2,3,4,7) 



Since Hi are written in terms of free fermions and x", it 
is straightforward to compute the right-hand side. The only 
non-zero OPE coefficients are listed below: 



Moreover, it is not difficult to verify that the above choice 
works for the equations for S2 as well: 

sin^ (0 r^' + cos^ (0 X^^x' + l^-^nei^x'-exn • 

(D3) 

Therefore, we may conclude that the following color- and chi- 
rality (L/R) independent SO(2) transformation R{9) generates 
the 'duality' in the continuum limit: 



9 9 
GR = CL/RCOs--XuRsin- , 

9 9 

XL/R = ?L/RSin-+xJ^/RCOS- . 



(D4) 



This is not the end of the story. Besides the uniform part 
discussed above, the continuum expression of the SU(4) gen- 
erators contain two kinds of oscillating parts; one has momen- 
tum 7r(= 4A;p) and oscillates in a staggered manner, and the 
others exhibits ±7r/2(= 2k-p) oscillation. The argument goes 
similarly for the 4fcF-part and we can verify that the choice 
(D4) works. 

The easiest way to write down the transformation of the 
2fcF-part would be to use the Dirac fermions instead of Majo- 
rana fermions ^ and x and look for the correct transformation 
rule for the Dirac quartet. Again, as in Eq. (Dl), we require 
that the 2fcF-part of the transformed SU(4) generators 

G = DX^R (2fcF) or R^X^L {-2kf) 

be correctly reproduced by linear combinations of the original 
generators as in Eqs. (9a) and (9b). The answer is quite simple 
and is given as follows: 



where 



R = U{9)R , L = U{9)L 



R = (*R,1T, *R,U, *R,2T, *R,2J,) 
L= («'L4T,*L4i,^L,2T,^'L,2i) 



(D5) 
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APPENDIX F: GROUND-STATE DEGENERACY 

In this appendix, we describe how to obtain ground-state 
degeneracy by inspecting the classical ground states within 
the bosonization approach. As has been mentioned before, 
bosonic expressions of the physical operators have a kind of 
'gauge redundancy' and the ground states should be counted 
modulo the gauge redundancy. 

To this end, we have first to identify physical operators for 
which we define 'gauge transformations'. One obvious choice 
may be U(3) Dirac fermions in Eq. (82). However, the expan- 
sion of physical (lattice) operators in terms of the continuum 
ones contains the 2fcF-terms which cannot be expressed by 
these Dirac fermions. This obscures how to impose the gauge- 
equivalence for three bosons (p, (f^, and (pf. Instead, we start 
from the four Dirac fermions used in Appendix C to obtain 
the effective Hamiltonian of the SU(4)-model. Throughout 
this section, we shall use the same notations as in Appendix 
C. By using the SU(4) bosons (l>s, $f, and $sf (see Eq. (C6)), 



the interactions Va 
read: 



-A^ {Oiy {A =Q, SD, SC, and RQ) 



2A 

'^'o = ^{ cos("s/47r$s) cos("s/47r$f) 



and U (9) is given by Eq. (10). 



cos(V47r$s) cos(v47r$sf) 
-I- cos(V47r$f) cos(V47r$sf) 



(Fla) 
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VSD 



2A; 



SD 



cos(v47r$s) cos(v47r$f 



+ cos(v 47r$s) cos( V 47r9sf) 



Vrq = - 



— cos( V 47r$f ) cos(v 47r0sf) 

(Fib) 

^ sin + e, + $f - Of)] 

X sin ($s - + $f + Of)] 
(T^ sin [Vtt ($s + B,s + $f - Of)] sin( ViTrOsf) 

sin([V^($s - e, + $f + Of)] sin(%/4^esf)| 

(Flc) 

^ ( - sin ($, + e, + <i>f - Of)] 

X sin [V^ - Bs + $f + Of)] 
a"- sin( [Vtt ($s + 65 + $f - 9f )] sin(\/47r$sf) 

cr"" sin [Vtt ($s - 6s + $f + 6f)] sin(V47r$sf)| ■ 

(Fid) 



These expressions can be readily obtained from Eqs. (83). 
The 2x2 Hamiltonian is diagonal and we can freely choose 
one of the eigenvalues (say, +1) of the Pauli matrix cr^ 
throughout the calculation. In principle, our system may be 
described only by these three bosons. However, the relation- 
ship between three bosons <l>s(6,s), $f(8f), <i>.sf(0sf) and phys- 
ical operators is not so obvious and it is convenient to go back 
to the original (2-band) Hubbard model (see Appendix C). For 
this reason, we recover the charge boson $c(0c) and add the 
following umklapp term to find the semiclassical vacua of our 



problem: 



-A, 



umklapp COS I 



(F2) 



The procedure is as follows. First of all, the gauge 
redundancy'^ of the original four Dirac fermions (see Ap- 
pendix C for the definition of these fermions) 



exp -iV47r<I>ax 



■\/27rao 

*a,R = — ^^CXp (+iV4^$a.R 

(a=l,...,4) 



reads: 



$a,L/R ^ ^'qX/R + \/^ A^, 



aX/R 



(A^ax/R e Z) . (F3) 



That is, physical quantities should be unchanged even if we 
make the above shift. Therefore, the semiclassical ground 
states which are 'gauge-equivalent' by this shift should be 
treated as one and the true ground states are the equivalence 
classes of this gauge transformation. 

For the two phases Q and SD, a straightforward semiclas- 
sical analysis (for either $ or 0) is applicable and we can 
proceed in essentially the same manner as in Ref. 15. For the 
cases of SC and RQ phases, however, the situation is slightly 
tricky. Since both $ and 8 appear in a single sine/cosine in- 
teractions, we should introduce new fields 



(F4) 



before applying a semiclassical argument. Then we can read- 
ily find the semiclassical ground states of Vsc.rq to apply the 
method of Ref. 15. 
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